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Abstract 

We compute the coefficient function for the charge-averaged iy*^-exchange structure function F3 
in deep-inelastic scattering (DIS) to the third order in massless perturbative QCD. Our new three- 
loop contribution to this quantity forms, at not too small values of the Bjorken variable x, the 
dominant part of the next-to-next- to-next-to-leading order corrections. It thus facilitates improved 
determinations of the strong coupling ttg and of l/Q^ power corrections from scaling violations 
measured in neutrino-nucleon DIS. The expansion of Fj, in powers of ttg is stable at all values of 
X relevant to measurements at high scales Q^. At small x the third-order coefficient function is 
dominated by diagrams with the colour structure d"'"^dahc not present at lower orders. At large x 
the coefficient function for F3 is identical to that of Fi up to terms vanishing for x — > 1. 



1 Introduction 



Structure functions in inclusive deep-inelastic lepton-nucleon scattering (DIS) are among the most 
important observables probing QCD, the theory of the strong interaction. Present data on these 
quantities can provide accurate information about the quark and gluon momentum distributions 
in the proton down momentum fractions x ~ 10^^, see Ref. [1]. These distributions, in turn, are 
indispensable ingredients for the analysis of all hard (high scale/mass) scattering processes at 
proton-(anti-)proton colliders, cf. Refs. [2,3], such as Tevatron and the LHC which will form 
the high-energy frontier of particle physics for at least the next fifteen years. Structure functions 
are also very well suited for precision determinations of the strong coupling constant Ug, one of 
the fundamental parameters of our description of nature, in the framework of perturbative QCD. 

Inclusive quantities such as structure functions are the observables best accessible to field- 
theoretic calculations to 'high orders' (today: two loops and beyond) in the expansion in powers 
of the coupling constant. Indeed, the second-order partonic cross sections (coefficient functions) 
for inclusive DIS were completed as early as 1991/2 [4-6], while corresponding quantities for jet 
shapes in e+e^ collisions have been presented only very recently [7-9]. At large x (with hindsight: 
X > 10^^) the former quantities are the dominant part of the next-to-next-to-leading order (NNLO) 
contributions in renormalization- group improved perturbation theory. They are not sufficient, 
though, for NNLO analyses over the full x-range opened up by HERA and required for the LHC. 

The three-loop corrections for inclusive DIS - required to complete the NNLO framework 
and to provide the dominant next-to-next-to-next-to-leading order (N^LO) corrections at large-x 
enabling a perturbative accuracy of 1% for ttg determinations - have been the subject of a long- 
term research program, which started from sum rules [10, 1 1], and proceeded via low integer Mellin 
moments [12-14] to the computation of the exact expressions for all NNLO splitting functions 
[15, 16] and the third-order coefficient functions for the structure functions Fi and F2 [17, 18]. 
In the present article, we extend the latter results to the most important structure function not 
covered by Ref. [18], the vector- axial- vector interference structure function F3 in charged-current, 
specifically (W^ + W^)-exchange DIS measured with high precision in neutrino-nucleon DIS [1]. 

The remainder of this article is organized as follows: In Section 2 we briefly recall the general 
formalism for the calculation of the coefficient function for F3 and discuss some aspects specific to 
the present computation. We then write down the coefficient functions for in Section 3. 

The second- and third-order quantities are presented via compact and accurate parametrizations. 
We also address the behaviour of the third-order coefficient functions at large and small x, stressing 
the importance of the d^^'^dabc contribution not present at lower orders. This and other numerical 
effects are then illustrated in Section 4 where we assess the size of the higher-order corrections. 
Finally we summarize our results in Section 5 and close with a brief outlook on possible future 
improvements on the present accuracy. Appendix A contains the very lengthy exact expression 
of our new third-order coefficient function, and Appendix B complements the discussion of F3 in 
Section 3 by providing the subleading large-x logarithms for the case of Fj. It turns out that the 
quark coefficient functions for F3 and 2xFi = F2— F^ in charged-current DIS are the same in the 
large-.x: limit, specifically that C3 (x) = Ci{x) + 0{l —x) holds to (at least) order a^^. 
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2 Formalism and calculation 



We are interested in unpolarized charged-current deep-inelastic scattering (DIS), i.e., the reaction 

l{k) + nucl(;?) ^ l'{k') + X , (2.1) 

where /, /' and 'nucl' denote a charged lepton, its (anti-) neutrino (in this or the opposite order) 
and a nucleon with respective momenta k, k' and p. X stands for all hadronic states allowed by 
quantum number conservation. The inclusive cross section for the process (12.11) can be written as 
do ~ L'"^Wuv in terms of leptonic and hadronic tensors L^y and Wuv The former tensor is well 
documented in the literature for both pure electromagnetic and electroweak gauge-boson exchange, 
see, e.g., Ref. [1]. Thus there is no need to consider it here. 

Our focus is on the (spin- averaged) hadronic tensor W^v which can be expanded to define the 
structure functions Fj. 3 viz 



^,v{P^q) = ^|rf'ze'^-^(nucl,p|7;(z)7v(0)|nucl,p) 



Here q denotes the momentum transferred by the V7-boson, with = —q^ > 0, and represents 
the weak current. The Bjorken variable is defined as x = / {2p -q) with < x < 1. We do not 
consider the functions F2,l corresponding to the symmetric e^y and d/jy in this article, see Ref. [18]. 
Instead we address the structure function F3 associated to the totally antisymmetric tensor £^va|3 
which arises from the vector/axial- vector interference of the two V—A currents. Specifically we 
are interested here in the charged- averaged quantity which has been accurately measured 

in neutrino-nucleon DIS off (almost) isoscalar targets [1] (recall that F^^^^ is small in this case, 
being proportional to a flavour asymmetry in the quark sea). 

As in the previous calculations in Refs. [15-20], we derive analytic expressions for the Mellin 
moments of the perturbative contribution to structure functions, which for the present case are 
defined as 

fW^+w-^j^^q2^ = (\xx''~'F^^+'^-{x,Q^). (2.3) 

J 

To this end we make use of the optical theorem relating the hadronic tensor in Eq. (12.21) to the 
imaginary part of the forward Compton amplitude T^v of virtual gauge-boson-nucleon scattering. 
This amplitude is expressed in terms of a time-ordered product of two local currents to which 
standard perturbation theory applies, 

T^y{p,q) = i j dW-^nuc\, p\t(^jJ,{z)MO)') \nuc\, p) . (2.4) 
In D = 4 — 2£ dimensions the structure function F^ is obtained from T^y using the projection 



;i-2£)(l-£) Ip-q 



P^^ = _i___l__eW--M _ (2.5) 
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As discussed in detail in the literature, see, e.g., Refs. [19,20], the operator-product expansion 



(OPE) can be applied to the product of currents in Eq. (12.41) together with a dispersion relation. For 
the structure function considered here this procedure finally yields 
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^3 



^3 (^'fl-'^s )Av,nucl(A^,i"^ 



(2.6) 



One thus obtains the odd-A'^ moments (12.31) which uniquely fix all moments, and hence the complete 
x-dependence, by analytic continuation. The notation for the coefficient function of '^^ 
reflects this basis of odd moments. Ay.nucKA^jA'^) in Eq. (12.61) is the nucleon matrix elements of the 
combination Oy of local spin-A'^ twist-two quark operators corresponding to the total valence quark 
distribution, renormalized at the scale fx. l/Q^ corrections have been disregarded in Eq. (12.61) . 

As usual, we perform the renormalization in the modified [19] minimal subtraction scheme 
[21] in dimensional regularization. Hence also the renormalized coupling in Eq. (l2.6l) refers to 
D = 4 — 2£ dimensions, i.e., its scale dependence is given by 



d 



J as 



£ as - Po - Pi al - P2 at - P3 a. 



(2.7) 



where p„ denote the coefficients of the usual four-dimensional MS beta function of QCD. The 
renormalization to the third order requires the coefficients up to P2 [22,23], and the consistent 
application of the resulting three-loop (N-^LO) coefficient functions in data analyses even P3, the 
highest contribution to Eq. (12.71) computed so far [24,25]. 

There is one delicate issue entering the operator renormalization which deserves special atten- 
tion. This is the presence of 75 due to the axial-vector coupling of the -boson. As before in 
Refs. [11, 14, 15,20] we employ the so-called Larin scheme [11,26] which, at the present level, 
is equivalent to, but computationally more convenient than the original prescription [27, 28] for 
consistently dealing with 75 in dimensional regularization. See also the discussion in Ref. [6]. 
Therefore we need to perform a special renormalization, based on relating vector and axial- vector 
currents, in order to restore the axial Ward-identities, 



Za Ov • 



(2.8) 



The constant Za for the axial renormalization and the finite renormalization Z5 due to this treatment 
of Y5 in the MS-scheme are known to three loops [1 1, 26]. The explicit expressions read 



Za 



Z5 



1+a^- 



1 

+ - 
£ 



22 4 

— — - CfHj- 



1 



484 2 176 16 2 

■ CfCA + -zrz- CpCi tif — — Cf-n 



27 



308 



3578 



— CpCA + ~gj~ Cp^A + -g CpUf — CpCArif -I- — CfUj- j 



832 



— ^f^A'if - — y^F'\f 
8 



(2.9) 



1 +as 



- 4C^ - 10£C^ + £2c^(-22 + 2 ^2) 



+ a" 



+ £ (132 -48^3) +C^Q 



7229 



22 



331 



107 



3 



+4 



c3(-^ + 96C3)+C^Q(^-160C3)+QCi(-^+56C3 



/ 5834 



2147 



2 / 62 32 X ^356 32 X 52 2 



(2.10) 



Both these expansions are given in terms of the renormalized coupling (12.71) . and Z5 is expressed to 
exactly the depth in e implemented in the present calculation (see below Eq. (12.211) for a comment 
on the positive powers of £ absent in Ref. [11]). 

We are now ready to briefly address the standard part of the renormalization of the spin-A'^ 
operators entering Eq. (12.61 ). 

Of" = ZvO,^^^ (2.11) 



where (as in some equations below) we have suppressed the dependence on N, and exploited the 
fact that Ov does not mix with other operators under renomalization. The corresponding (non- 
singlet) anomalous dimensions Yv, governing the scale dependence of the operators 0^, via 



JlnA/2^^ " ^^^^ ' 
are connected to these renormalization constants in the standard way, 



Yv 



dln/u- 



■ Zy j Zy 



1=0 



(2.12) 



(2.13) 



,(0 



The reader is referred to Ref. [15] for the relation of the expansion coefficients Yv to those of 
other non- singlet combinations. In terms of these quantities, the MS renormalization factor Zy in 
Eq. (|2.11|) is given by the Laurent series 



1 + as-y^y^ + al 



1 



,(0) 



1 



,(1) 



^{(yI«^-2Po)(yI«^-Po)yI«^ 
^{3YlV)-2PoYi^^-2PiY? 



2£ 



+ 3e^^ 



(2.14) 



Of course, the perturbative calculation of the anomalous dimension (12.121) and the coefficient 
function cannot proceed via Eq. (12.61) including the non-perturbative nucleon states |nucl,p). 
However, as the OPE represents an operator relation, the calculation can be performed using quark 
states instead. Hence we apply the Lorentz projector (12.51) to the forward virtual-M^-quark Comp- 
ton amplitude T^n analogous to Eq. (12.41) and obtain the perturbative quantities 13 which can be 
expanded in powers of the renormalized coupling = a^/ (47l) at the scale /j^ = Q^, 

T,iN) = {rf\N)+a,Tl'\N)+a^,rl^\N)+a',Tl'\N) + ..)jAy, (2.15) 

where the A''-independent Ay denotes the quark operator matrix element. After performing the 
operator renormalization according to Eqs. (12.91 ) - (I2.11|) and (|2.14l) . one thus arrives at explicit 
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expressions for the T^^^ . The anomalous dimensions yl,'^ can then be read off from the poles in 
£, while the coefficient function C3 is related to the finite terms in £. During this procedure it is 
important know the expansion of the D-dimensional coefficient function 

00 

C3-(as,A^) = Y^a'J4'\N)+e4'\N)+e^4'\N) + ...) (2.16) 

to a sufficient depth (positive powers) in £ at lower orders. We normalize the contribution to 
Eq. (12.151 ) including Ay, i.e., 

ll^iN) = 1 (2.17) 

implies 

cfliN) = 1 , 4'^(A^) = ^.('^^(A^) = ... = 0. (2.18) 

In this normalization also the first-order expression is the same for and ^3^"-*^-. 

Again suppressing the A/^-dependence it reads 

^il' = jtS' + 4!i + e4!1 + '='41 + o{e') . (2.19) 

The corresponding second- and third-order results are given by 

.(2) . 



+ 4'l+4:i4' + e{<l+4!lTS'} (2.20) 



and 

^(3) 



^{(tS'-2p„)(,S'-p„),S'} 

+ ^ {3yi" yS' - 2Poy<" - 2P, tS' + 34'i (tS' - Po) tS' 
+ i {2T?' + 34:iy<V64'lyS» + 34'i (yS>-P„)TS'} 

+ ct 44!iTt" +4'^yS> + i4'i (tS' - Po) yS' . (2.21) 

Consequently the two- and three-loop coefficient functions can be read off from the £-independent 

(2) (3) 

parts of T-^ _ and T-^J^ after subtracting the respective contributions due to the lower-order £ and 
£^ quantities in Eq. (12.161) . 

We are now in a position to give the comment announced above on the £-terms in Eq. (12.101) . 

(2) (2) 

These terms do not affect the extracted coefficient functions C3 _ and C3 _. Without them, however, 

the functions a^2±' ^3^1 ^^'^ ^^3- would exhibit an unphysical behaviour of their A/^-independent 
terms (which should be the same as those for F2), a feature irrelevant here but unwanted for more 
general applications such as the determination of time-like (fragmentation) coefficient functions 
via a suitable analytic continuation of the DIS results, cf. Refs. [29,30]. In other words, Eq. (|2.10l) 
effectively restores the anticommutativity of 75 also for those £ and £^ contributions. 
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The actual computation of the Feynman diagrams for the contributions to Eq. (12.151) follows 
those of Refs. [15-18,31] in every respect, so we can be very brief here. The graphs have been gen- 
erated automatically with the diagram generator Qgraf [32]. All further symbolic manipulations 
have been performed in FORM [33,34], using the Summer package [35] for the analytic evaluation 
of all nested sums. The calculation relied on a massive tabulation of intermediate integrals, and 
check of all intermediate and final results were performed for fixed values of against the Mincer 
program [36,37] and the findings to Ref. [14]. The three-loop 1/e terms of the present calculation 
have already been used to complete the set of non-singlet NNLO splitting functions [15]. 

There are two aspects which deserve special attention in the context of the present calculation. 
The first is the appearance of two functions g^, g2 in the final odd-A'^ results for the coefficient 
functions (see Refs. [18] for very similar even-A'' functions in the results of F2 and Fi) which fall 
outside the class of simple harmonic sums [35] sufficient at previous orders, 

gn{N) = A^" (5^5 -25_5 + 45-2^3 -45-2,-3 + 85-2,-2,1 +453,-2-454,1 +255) • (2.22) 

Note that the bracketed combination of harmonic sums vanishes as 1 /N^ for N ^ 0°, hence gn{N) 
do not contribute to the leading In'^A'^, = 1, . . . , 6, behaviour of the coefficient function c^^^_(N). 
The x-space expressions corresponding to Eqs. (12.221) can be found at the end of Appendix A. 

The second interesting point is the presence of a new colour structure, a contribution propor- 
tional to the higher SU(nc) group invariant d"'^'^dabc in the three-loop coefficient function c^^K This 
new colour factor is related to Feynman diagrams of a particular topology concerning the fermion 
flow through the diagram. We distinguish the two cases (so-called flavour classes) fl2 and //02 
depending on whether both VF^ -bosons are attached to the open fermion line of the initial and final 
state quark (fh) or whether both gauge bosons are coupled to a closed fermion loop {fhi). Sample 
diagrams for the two cases are displayed in Fig.[T] The impact of the new //02 contributions with 
d"-^'^dabc on the coefficient function will be discussed in the next sections. 




+ 348 other diagrams +11 other diagrams 

Figure 1 : Typical Feynman diagrams for the contributions without (left) and with (right) d^^'^dabc 
to the charged-current DIS coefficient function C3 up to the third order in the strong coupling a^. 
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3 Results and discussion 



We are now ready to present the charged-current coefficient function C3 to the third order in the 
strong coupling = a^l [An), i.e., the coefficients ci'f in 

Here ^vai represents the total (flavour summed) valence quark distribution of the hadron, ^vai = 
LiLiWi '~ Qi) where n^- is the number of effectively massless flavours. ® denotes the Mellin con- 
volution which turns into a simple multiplication in A^-space. All results below will be given in the 
MS scheme for the standard choice fi^ = /uf = of the renormalization and factorization scales. 
The complete expressions for the dependence on /j^ and /j^ to the third order can be found, for 
example, in Eqs. (2.16) - (2.18) of Ref. [38]. 

As discussed above, our calculation via the optical theorem and a dispersion relation directly 
determines the coefficient function (13.11) for all odd-integer moments N in terms of harmonic 
sums [35,39,40]. From these functions the jc-space expressions can be reconstructed algebraically 
[20,43] in terms of harmonic polylogarithms [41-43]. As in the case of F2,l in electromagnetic 
DIS presented before [17, 18], the exact third-order expressions are unpleasantly long in both 
A^-space and Jc-space. We therefore refrain from writing down the former in this article, and defer 
the latter to Appendix A. 

For the convenience of the reader we first recall the known results up to the second order. The 
first-order contribution to Eq. (|3.1I) is identical to that for the W^ — W^ case C^^ and given by [19] 

4%{x) = C^{42)i-3©o-(9 + 4C2)5(xi)-2(l+x)(Li-Lo) 

-4x-iLo + 4 + 2x} (3.2) 

with Cp = {n^ — 1 ) / (2nc) =4/3 and Hc = = 3 in QCD. Here and below we use the abbreviations 

Xi = \ —X , Lq = \nx , Li = Inxi , ©yt = [x^^Li]+ . (3.3) 

As usual, the -i-distributions are defined via 

f'dxa{x)+f{x) = Cdxa{x){f{x)-f{\)} (3.4) 

JQ Jo 

for regular functions f{x). Convolutions with the distributions 'D^ in Eq. (13.31) can be written as 

xl^..«/lW = -./w} l^^^^il^ . (3.5) 

Already at two loops the coefficient functions involve polylogarithms and Nielsen functions, 
hence it is convenient to have at one's disposal accurate parametrizations in terms of the elementary 
functions in Eq. (13.31) . After inserting the QCD values of the colour factors Cp and C^, the exact 
two-loop coefficient function (A. 3) first obtained in Refs. [6,20] can be represented by 
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c\ '_{^) ^ 128/9 ©3 -184/3 ©2 -31.1052©! + 188.641 ©o-338.5725(jCi) 
-16.40L^ + 78.46L?-470.6Li- 149.75 -693.2x + 0.218xL^ 
+ LoLi(33.62Lo- 117. 8Li)-49.30Lo-94/3L^- 104/27 

+ Hf { 16/9 ©2 - 232/27 © i + 6.34888 ©o + 46.8464 5(xi) + 0.066 
-0.663L? + 24.86Li -5.738-5. 845x-10.235x2-0.190xL^ 
+ 4.265 L0L1+ 20/9 Lo(4 + Lo)} (3.6) 

with an error well below 0.1% at all values of x. This expression is slightly less compact, but 
considerably more accurate than the (practically sufficient) previous parametrization in Ref. [44]. 
At this order the coefficient functions for the V7++W^ and VF+— W^^cases are different. The two- 
loop coefficient function c^'j^ for the latter can be evaluated via Eq. (13.61) above and Eq. (2.9) in 

Ref. [45] for 5c = cf^+-cfl_ . 

The coefficients of the +-distributions 'D^ in Eq. (13.61) and Eq. (13.71) below are exact up to a 
truncation of the values of the Riemann ^-function. Also exact are those coefficients of Lq = In^x 
and Lj = ln'^(l —x) given as fractions. Most of the remaining coefficients have been obtained by 
fits to the exact coefficient functions at 10^^ < x < 1 — 10^^. Finally the coefficients of 5(1 — x) 
have been adjusted very slightly from their exact values using the lowest integer moments, thus 
fine-tuning the convolution with the quark distributions to maximal accuracy (cf. the discussion at 
the end of Section 4 of Ref. [16]. 

We now turn to our new three-loop results. Inserting, as in Eq. (13.61) above, the numerical QCD 
values of the n^^-independent SU(nc) colour factors, the third-order contribution (A.4) to Eq. (13.11) 
can be approximated by 

c^^l(jc) ^ 512/27 ©5-5440/27 ©4 + 501.099^3 + 1171. 54©2-V328.45©i 
+ 4442.76 ©0- 9172.68 5(jci)- 512/27 Lf + 8896/27 1396 
+ 3990Lf + 14363 Li- 1853 -5709x + jcxi(5600-1432x) 

- LoLi (4007 + 13 12 Lo) - 0.463 ;c4- 293.3 Lo- 1488 L^- 496.95 4 

- 4036/8 IL^- 536/405 

+ {640/81 ©4 - 6592/81 ©3 + 220.573 ©2 + 294.906 ©1 - 729.359 ©0 
+ 2575.46 5(;ci ) - 640/8 1 ij + 32576/243 L\ - 660.7 + 959. 1 Li 
+ 516.1 -465.2;c + ;c;ci(635.3 + 310.4x) + 31.95xiL| 
+ LoLi(1496 + 270.1Lo-1191Li)-1.200xL^ + 366.9Lo + 305.32L^ 
+ 485 12/729 + 304/8 1 } 

+ {64/81 ©3-464/81 ©2 + 7.67505 ©1 + 1.00830 ©q- 103.2602 5(;ci) 

- 64/8 1 + 992/8 1 - 49.65 Li + 1 1 .32 - xxi (44.52 + 1 1 .05 x) 

+ 51.94X + 0.0641 xL^-LqLi (39.99 + 5.103Lo- 16.30Li) - 16.00 Lq 

- 2848/243 - 368/243 } 
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+ //o2 nf {2. 147 - 24.57 Li + 48.79 - {242 A - 150.7 x) - LqLi (8 1 .70 
+ 9.412Li) +xLo (218.1 + 82.27 L^) - 477.0Lo - 1 13.4L^ + 17.26 



16/27L^};c 



1 ■ 



(3.7) 



Here the factor //02 (= 1 for the numerical evaluation) indicates the d^^^dabc contribution entering 
at this order for the first time, cf. Ref. [11]. Also Eq. (13.71) . first presented in Ref. [46], deviates by 
much less than one part in a thousand from the corresponding exact expression which we evaluated 
using a weight-five extension of the Fortran package [47] for the harmonic poly logarithms. 
Eqs. (13.61) and (13.71) can be readily transformed to Mellin space at complex values of A'^ (using, e.g., 
the appendix of Ref. [40] for the moments of ln;c ln^(l —x) etc) for use with A^-space programs, 
such as QCD-Pegasus [48], for the evolution of parton densities and structure functions. 



The x^ \ and x^O end-point behaviour of the higher-order contributions to Eq. (13.11) is of 
special interest, both theoretically and phenomenologically. We first address the large-x limit. Here 
the leading terms of C3 j^{x) are the soft-gluon -i-distributions Ti^ with = 0, . . . , 2n — 1. To order 
ttg the corresponding coefficients are identical to those for the electromagnetic and charged-current 
structure functions F2 which can be found in Eqs. (4.14) - (4.19) of Ref. [18]Q 

The terms with 5(1 —x) (arising from virtual corrections and soft-gluon contributions) are the 

(3) 

same for all four charged-current coefficient functions 3 _|_ and differ from the corresponding 
photon-exchange quantity of Ref. [18] only by the obvious absence of contributions from the fl\\ 
flavour classes (where the photons couple to different quark loops) in VK-exchange processes. For 
the convenience of the reader we here collect the third-order 5(1 — x) terms scattered over the 
13 pages of Eq. (A.4): 



.(3) 
'3,± 



6(xi) 



C^Cp 



1909753 



1944 



143255 ^ 105712 ^ 25184 ^2 416 „ 

C2 + ^:^C3 + -r^C2--^C5 



+ 540^2^3 

6419 



248 3 



r 87632 ^2 

C3 + ^C2 

7255 3379 



81 
3512 2 

4952 



81 

+ QC^ 
6644 



135 

9161 104117^ 



9 

2148 



12 54 
^1016 3 33556 2 

C2C3 + ^C2-3Yr^3 



24 



^2 - 3 1 8 C3 - ^ C2 + 1 240 ^5 + 808 C2C3 



304 3 4184 2 
'=2 + oi c 



315 



+ Cpn 



+ Cj^Cpnf 



+ CpUf 



142883 40862 ^ 



5 

_ 9517 
~l86" 
18314 



860 

~2n 



152 



^3-9^1 



486 
341 



81 



81 



2488 84, ^ 4, 

C2 + tC5-^C2C3 



135 



36 



5491 ^ 1348 . 

^C2 + ^C3 



16472 2 592 
135 ^2~~9~ 



3 



32 
27 
56 

352 



C2C3 



(3.8) 



^ The sign of the term —232 ^5 in the first line of Eq. (4. 1 9) has been misprinted in the hep-ph version of Ref. [18]. 
The journal version is correct. 
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Finally we consider the subleading (integrable) large-x logarithms which have attracted re- 
newed theoretical interest recently [49-51]. Terms up to ln^"^^(l —x) occur in the n-th order 
coefficient functions. Their coefficients are the same for the W+-f and —W^cases. Unlike 
the -i-distributions, however, these contributions differ (except for the highest power in each colour 
factor, where the coefficient is minus that of the highest -i-distributions) between F2 and F3. For 
the present structure function the two-loop contributions read 



.(2) 
■3,± 

.(2) 
■3,± 



.(2) 
■3,± 



22 



Cj^Cp + 52 C 



640 



8C2 



124 



C/[16-32C2] + -^CpUf 



The corresponding coefficients for the third-order coefficient functions are given by 



.(3) 
■3,± 

.(3) 
■3,± 

.(3) 
■3,± 



8C. 



220 



40 



+ 84C;? ~ —Cpn 



484 
176 



cic 



9056 
27 



.(3) 
■3,± 



27 



C^Cp 



1640 



■CpHf 



32^2 
16 



C3[110-96C2] 



27 



pnf 



7580 98 
27 ^T^^ 



^F 



1097 



+ 656^2 + 16^3 



- CpHf 



2098 112 



+ 



248 
27" 



12031 
9 

/-^ 2 
■CpHj 



-372^2-240^3 
2734 16 „ 



27 



.(3) 
■3,± 



— C^Cp 
_ 12568 

+480^3 

+ Cjtljr 



138598 4408 ^ ^^^^ 176^2 
-^C2 -272^3 + ^^2 



81 



1904 ^764 2 
— ^3 + 4,2 



376 



+ C'p 
45260 



1741 



+ ■ 



Cj^Cp 
1220 



69833 
162 



81 



6 3 
108^2- 16^3 



C2 



9763 2224 



81 



^2 



112 



^3 



Cptij- 



3520 _ 32 ■ 

81 "y^^ 



(3.9) 
(3.10) 
(3.11) 

(3.12) 
(3.13) 



(3.14) 



(3.15) 



(3.16) 



The extraction of some of these coefficients from Eqs. (A. 3) and (A.4) is far from trivial. We there- 
fore provide the corresponding results for F2, which we did not include in Ref. [18], in Appendix B 
where we also discuss an unexpected (to us) relation between the large-x coefficient functions. 
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We now turn to the small-x limit of the second- and third-order coefficient function in Eq. (13.11) . 
The leading terms in this case are the 'double-logarithms' Lq = In^x with k — 1, ...,2n — l at the 
n-th order in ttg. The two-loop coefficients are 



.(2) 

'3,- 



.(2) 

■3,- 



.(2) 

■3,- 



2 C jiCp -\- — Cp 



1 
3 



4 



^-0 



103 2 5 

Cj^Cp + 2\Cp + -CpHf 



6 

Cji^Cp 



122 



8^2 



20 



+ [21 + st;2] + —CpHf . 



(3.17) 
(3.18) 
(3.19) 



(2) 

The results for the corresponding W^—W' coefficient function ^ can be obtained by combining 
Eqs. (I3J71) - (I3J91) with Eq. (2.7) in Ref. [45]. Also that equation for 5c 



(2) 



.(2) 



C3 i contams 



terms up to In hence c\ _ and Cj^^ already differ in the leading logarithm. On the other hand, 
the leading logarithms (but only these) are the same for the e\en-N based two-loop coefficient 

(2) (2) 

functions C2 ^ and C3 ^ for F2 and F3. 

At the third order in ttg, the coefficients of In^ jc in Eq. (13.11) are given by 



.(2) 



.(3) 
•3- 

.(3) 

■3,- 



29 



53 



32 d^'^'^dabc 



+ 5^^^^ " 15^-^^^ + 30^^ " 15 n, 

166 2^ 43 2 89 3 46 
^ + ^ ~^ T2 27 



31 

Is 



■CpHf 



(3.20) 
(3.21) 



.(3) 



.(3) 

'3 - 



9799 64 „ 



-Ci 



49 710 .' 

T + ^^2 



+ 



772 

-^C^CpUf 



24245 220 „ ' 



2179 



■CpHf 



92„ 2 
81^^"/ 



704 256 



— C^Cp 



9 3 
38642 226 



^2 



81 
13297 



r 212^ 
^^2 + — ^3 



4216 112 



27 



54 

1085 1373 



6 

^ ^abc 



+ 401^2 + 328^3 
^2- 286^3 



Cpfif 



2687 62 „ 



712 



1056-96^2-^^3 



(3.22) 



(3.23) 



.(3) 

■3,- 



53650 11260 ^ 1696 

^2 + ^ 



81 



27 



416 

y 



C3 + ^Cf 



C^Cp 



404119 
324 



11 



16273 



27 
1978 



^2 



Cpflj: 



^3- 

"271 



9206 
4748 



C3 



3416 2 



5553 2308 



^2 



15 
718 



+ C^Cf 



26084 
81 



J62 



+ 



27 



1028 ^ 

^2 + ^^C3 



^2 



3112^ 304 / 
1684 



f '7 



81 



16 . 
yC2 



1600 



1024 



^2 



1088 



^ 2176^ 

^3 + ^^: 



(3.24) 



Inserting Ca = 3, Cf = 4/3, d'^^'^ date / = 5nj/l% and the numerical values of the ^-function, 
Eqs. (13:201) - (13:241) yield 



.(3) 
■3,- 

.(3) 

■3,- 

.(3) 
■3,- 

.(3) 

■3,- 

.(3) 

■3,- 



1.32346- 0.59259 n^//02 



-49.8272 + 3.75309 n 
-496.842 + 66.5460 n 
-1485.12 + 305.353/1 
-274.294 + 367.303/1 



1. 5 1440 /i/+ 17.2626 /2^/Zo2 
114.125/2^/Zo2 
494.486 /i^//o2 



n.llQlnj 



160.229 /i/ 



(3.25) 



where, as in Eq. (13.61) above, //02 marks the d"^'^dabc contribution. 



The qualitative pattern in Eq. (13.251) is the same as found for the three-loop splitting functions 
in Eqs. (4.16) and (4.18) of Ref. [15] and the non-singlet photon-exchange coefficient function 
for F2 in Eq. (4.26) of Ref. [18]: The coefficients show a dramatic rise until the fourth term, 
thus precluding any meaningful approximation by (up to) the three highest contributions (13.201) - 
(13.221) at practically relevant values of x. It is also worthwhile to note that the new d"^'^dabc colour 
structure actually forms the largest contribution to the first line of Eq. (13.6!) for physical values of 
the number of flavours. Hence it is not possible to obtain, or reliably estimate, even the leading 
a" In^"^^ a: terms at higher orders n by any resummation of lower-order structures. 

For the corresponding — W^ coefficient function c^^^_^ only the six lowest even-integer 
moments have been computed so far [52,53]. As observed below Eq. (13.191 ), the leading-log 
coefficients are the same for the even-N coefficient functions for F2 and Ft, to two loops. Assuming 
that this relations holds at three loops as well, we can employ Eq. (|3.20l) above and Eq. (4.14) in 
Ref. [18] to derive the conjecture 



5c 



(3) 



1 



(Q-2C/c2--(Q-2C/2c^ 



F ■ 



(3.26) 



The relation is supported by the fact that, as the even-integer moments [45] of 5c 



(3) 



.(3) 

■3,-F ■ 



.(3) 



(with the Cj^-specific d"^^'dabc contribution removed from the latter quantity), also the coefficient 
(|3.26i) shows the characteristic l/n^ suppression in the limit of a large number of colours Hc pre- 
dicted in Ref. [54]. We expect to be able to verify the above conjecture by a complete calculation 

(3) 

of the coefficient function c\ {{x) in the not too distant future. 
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4 Numerical implications 



We start this section by graphically illustrating our new third-order coefficient function c 3 _ (x) and 
its convolution with a schematic but sufficiently characteristic quark distribution, 



xf{x) 



,0.5/ 



1 



(4.1) 



All curves in the corresponding two figures are scaled down from the normalization (13.11) by a 
factor 2000 ~ (47r)^, effectively switching back to the normal-sized expansion parameter . 

On the left side of Fig. [2] we compare the exact x-shape given by Eqs. (13.71) and (A. 4) to 
two approximations indicating the previous uncertainty band [55]. This band was constrained by 
the lowest seven odd-integer moments A'^ = 1 ... 13 computed in Ref. [14] and the four leading 
-i-distributions [56] fixed by the next-to-leading logarithmic threshold resummation matched to 
the second-order coefficient function of Ref. [6]. We see that this estimate was indeed reliable, 
and sufficiently accurate atx^ 0.2. Thus a future determination of the fourth-order effects can be 
started reliably from fixed-// moments combined with the seven leading -i-distributions now fixed 
by the next-to-next-to-next-to-leading logarithmic threshold resummation performed in Ref. [57]. 
We refer the reader to this article for a discussion of higher-order effects at very large values of x. 
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Figure 2 : The third-order coefficient function c 3 l (x) for four flavours , multiplied by ( 1 —x) / 2000 
for display purposes. Also shown (left) are the previous uncertainty band [55] and (right) the 
corresponding contribution to the H^^+VK^ structure function F2.ns (from Ref. [18] for flu = 0) 
and the small-x approximations by the leading and next-to-leading logarithms (13.201) and (13.211) . 
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The right side of Fig.[2lfocuses on medium to small values of x. A huge low-x rise is found from 
X ~ 10^-^. This is totally different from the behaviour of the corresponding electromagnetic and 

+ coefficient function for F2 also shown in the figure. In fact, this rise can be attributed 
entirely to the new d"^'^dabc colour structure in Eqs. (13.71 ) and (13.251 ): If this contribution were 
removed from the coefficient function, then the solid line (red in the archive version) would not 
end at 36.4 for x = 10^^, but at —4.3, in qualitative agreement with the behaviour of C2^^ix). 

(3) 

What physically matters, of course, is not the distribution (aka generalized function) C3 l(x) 
itself but the resulting contribution to the structure function (13.11) . obtained by the convolution (13.51) 
and its obvious counterpart for the regular terms with the valence quark distribution. This convo- 
lution is shown in Fig.[3]for the typical quark distribution (14.11 ) (the normalization is irrelevant as 
we display all results normalized to /, thus suppressing large but trivial variations over the chosen 
wide range in x). For the sake of a direct comparison the same shape is used in Fig. [3] for the 
(different) quark distribution entering F2,ns • 

(3) (3) 

A comparison between the result for Cg _{x) and l^^{x) illustrates the 'x-shifting' power of 
the convolution integrals: The former coefficient function (for //02 = 0) is larger than the latter at 
all X < 0.8, yet its convolution result is larger only for x < 10~^. Note also that the small-x rise 
for F3, already delayed in Fig. [2]by about one order of magnitude in x to x < 10 by non-leading 
logarithms, is confined to x < 10^^ after the convolution. 
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Figure 3: Convolution of the V7"''-|-V7 charged-current coefficient functions c\ _ and for 
= 4 with a simple shape typical of non-singlet quark distributions of the nucleon. Also shown 

(dashed) is the result obtained for the former if the d"^^'dabc part of Eqs. (13.71) and (A. 4) is left out. 
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The normalized convolution {c\^ _® f) / f of Fig. 2, i.e., the coefficient of in +^ /^vai 
for x^vai — x^^^{\ —xY, exceeds three, thus 2.5% of the lowest order for tts — 0.2, only at x > 0.75 
(where we run out of measurements far from the resonance region anyway) and x < 1.5 ■ 10^^. 
Hence the perturbative expansion for this structure function appears stable at all practically relevant 
values of x, including the very low values x ^ 10^^ probed at ^ lO'^ GeV^ by the scattering of 
ultra-high energy cosmic neutrinos, see, e.g., Ref. [58]. 

We finally assemble Eqs. (13.21) . (13.61) and (13.71) to illustrate the perturbative expansion (13.11) of 
the structure function +^ . In order to easily compare the contributions of the various orders, 
we use an order-independent value of the strong coupling, a^ifjf = /jf = Q^) = 0.2, corresponding 
to a scale 2^ ~ 30. . .50GeV^. For the same reason we use the same quark distribution at all 
orders of Eq. (13.11) , starting with the model shape (14.11 ) already employed above. 

In Figs. |4] and [5] we show the total results at N'LO as well as the relative effects of the /-loop 
contributions to the structure function. The higher-order (NNLO and N-^LO) corrections are small 
(1% or less at N^LO) except for the soft-gluon rise towards x= I. Note that with increasing order 
this rise becomes steeper, but starts at larger values of x. For instance, the relative /-loop corrections 
exceed 5% atx > 0.46, 0.7, 0.83 for / = 1, 2, 3, respectively. The latter value corresponds to an 
invariant mass W GeV for = 40GeV^. 
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Figure 4: The structure function (13.11) to order (N^LO) for a fixed value of tts and the schematic 
valence quark distribution (14.11) . Next to the cumulative effect of the known corrections on the left, 
we show on the right the relative corrections due to the /:-loop coefficient functions, including an 
N'^LO large-jc estimate by the -i-distributions ^^i,...,? known from the threshold resummation [57]. 
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Figure 5: As Fig. HI but for the higher-order corrections for down to small values of x. 



An estimate of the large-x fourth-order coefficient function by the seven -i-distributions fixed 
the next-to-next-to-next-to-leading logarithmic soft-gluon (threshold) resummation [57] strongly 
suggests that this trend will continue at even higher orders. As shown in the right part of Fig. HI 
this estimate leads to x > 0.9 for a relative four-loop correction exceeding 5%. All in all, the N-^LO 
expansion can be considered safe to, at least, x ^ 0.8 at 40GeV^ (the safe range, of course, 
widens (shrinks) with increasing (decreasing) due to the scale dependence of tts). 

Specific numbers as given in the last two paragraphs depend on the quark distribution. This 
is illustrated in Fig. [6l where the Reggeon x counting-rule ansatz (14.11) . cf. Ref. [59], for x^val 
is modified at large x (left) or small x (right). Suppressing ^vai in either region leads to larger 
corrections mainly in the same region, hence Fig. [6] focuses on the large-x region in the former 
case, and small values of x in the latter. Again we show the results for a, = 0.2 and four flavours. 

Suppressing x^yai by two powers of ( 1 — jc) at large-jc leads to a considerable widening of the 
region of large soft-gluon corrections, with the NLO, NNLO and N^LO 5% x-values listed above 
Fig. m all reduced by about 0. 1 . A reduction of the small .)c-quark distribution by a factor x'^^^, i.e., 
a factor of 10 at jc = 10^^ with respect to Eq. (14.11 ) has a rather dramatic effect, especially at NLO, 
as shown in the right part of Fig.[6l However, while the relative NLO correction is as large as 24% 
dX X = 10^ in this case, the higher-order corrections remain small, with the N-^LO contribution 
exceeding 1% only at x < 2 ■ 10^^. Hence, even under these conditions, the perturbative expansion 
to N^LO proves sufficient for (more than) all practically relevant situations. 
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Figure 6: As the right parts of Figs. |4] and [51 but for different choices of the quark distribution ^vai • 

Like the large-x rise in the right part of the figure, the small-jc rise on the left appears to move 
closer to the end-point as the order in tts is increased. Unlike at large x, though, there is no way 
of estimating the fourth-order contributions at values x < 10^^, as there was no way predicting 
the size of the d"^'^dabc three-loop correction from lower-order information. Also fixed moments 
cannot provide constraints at such values of x, therefore the present (fortunately satisfactory) status 
will remain unless / until someone (else) calculates the exact fourth-order coefficient function. 

5 Summary and outlook 

We have extended our previous computations [17, 18] of exact third-order coefficient functions in 
inclusive deep-inelastic scattering to the charged-current structure function +^ . Hence the 
next-to-next-to-next-to-leading order coefficient functions are now known, in massless perturbative 
QCD, for all structure functions for which precision measurements have been performed in fixed- 
target DIS and/or at HERA, enabling improved analyses of such data at jc > 0.01. 

Also the present calculation has been performed in Mellin-A'' space, obtaining an analytic for- 
mula in terms of harmonic sums [35] for all odd-N moments (as in Refs. [19] and [20] at first and 
second order) using the optical theorem and a dispersion relation in the Bjorken variable x. From 
this result, which we will make available but did not write down in this article for brevity, we have 
reconstructed the equally lengthy exact jc-dependence in terms of harmonic polylogarithms [43] 
presented in Eq. (A. 3). A compact and accurate parametrization is provided by Eq. (13.71) . 
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The singular large-x terms, (1 — ^)_|_^ ln^(l —x) and 5(1 —x), are the same for and 
pW+±w rpj^g latter coefficient function, in turn, is identical to that for Fj® ™' of Ref. [18] up to 
the obvious absence of the flu diagram class which contributes to the coefficient of 5(1 —x) at 
three-loop in the photon- and Z-exchange cases. Even at the largest values of x practically relevant 
at large scales and invariant masses W^, the observable- specific non-leading contributions to 
the coefficient functions are non-negligible. We have presented, in particular for use in theoretical 
studies of subleading terms, explicit expressions of the ln'^(l —x) contributions to both F2 and F3 
which turn out to be related via the corresponding terms for already presented in Ref. [18]. 

At small-x the N'LO non-singlet coefficient functions include potentially large logarithmic 
terms up to In^'^'x. For Fj*^ ™' we found that the prefactors of the third-order terms are such 
that a small-x rise only occurs at irrelevantly low values of x [18]. The results for the vector- 
axial-vector interference quantity F^ '^^ would be similar, were it not for the d"^^dabc fhi 
diagram class (absent in F2 due to Furry 's theorem) which occurs at the third order for the first time. 
These diagrams dominate the small-x limit and lead to a rise of the coefficient function at jc < 10 
before and x < 10^^ after the convolution with the valence quark distribution. Nevertheless, the 
corrections remain very small at jc-values accessible to colliders and unproblematic down even to 
the very low values of x of interest in DIS of ultra- high energy cosmic neutrinos at 2^ ^ M^. 

Progress beyond our present three-loop accuracy for p^^^^ would be possible at large-x if the 
fixed-A'' results of Ref. [14] could be extended to the fourth order, by combining those results with 
the substantial constraints from the threshold resummation [57] into approximations analogous to 
those of Ref. [55]. Given that a first four- loop moment has been computed already, the splitting 
function for quark combinations such d& u + ii — {d + d) atN — 2 for three flavours [60] , this seems 
a not entirely unrealistic perspective. 

On the other hand, any attempt of reliably inferring the small-.x: behaviour of f^'^^^ via a 
resummation of lower-order information appears doomed to failure by the hierarchy of the small-x 
coefficients, recall Eq. (|3.25l) . and the possible occurrence of dominant new colour factors such as 
d'^^'^dabc in the present three-loop case, see also Ref. [15]. Hence the rather forbidding extension 
of the present all-A'^ computation to the fourth order would be required for progress at small x. 
Fortunately the size of the three-loop corrections does not call urgently for such a calculation. 

Form files of our results in both A'^-space and x-space, and Fortran subroutines of the exact 
and approximate coefficient functions can be obtained from the preprint server http : // arXiv . org 
by downloading the source of this article. Furthermore they are available from us upon request. 
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Appendix A: The exact x-space results 



Here we present the exact expressions for the coefficient functions up to the third order in terms 
of harmonic polylogarithms H,„j ....,„j^,(jc), mj = 0, ±1 [43]. Functions up to weight (number of 
indices) 2n — 1 contribute at order a". For 'brevity' we suppress the argument x and define 

PqqW = 2{l-x)-'-l-X. (A.l) 

All divergences for x I in Eq. (lA.ll) and below are to be read as +-distributions, see Eqs. (13.41) 
and (13.51 ) above. In this notation the one-loop coefficient function (13.21 ) for can be written as 

4%{x) = C^(^(5+x)-^77qqW(3 + 4Ho + 4Hi)-5(l-.^)(9 + 4C2)) . (A.2) 

The exact two-loop result corresponding to Eq. (13.61) reads 

4^l{x) = C^(-4pqq(-x)(7^3-8H_i^2-2Ho + 2HoC2-8H_i,_i,o + 10H_i,o,o 

+ 4H_ 1,0, 1 + 6Ho,- 1 ,0 - 3Ho,o,o - 2Ho,o, i ) + ^T^qq W (5 1 + 1 28^3 + 48^2 + 1 22Ho + 96H0C2 
+ 54Hi + 32H1C2 - 12Ho,o - 48Ho,i - 72Hi,o - 72Hi,i + 96Ho,_i,o - 32Ho,o,o - 96Ho,o,i 
- 96Ho,i,o - 1 12Ho,i,i - 48Hi,o,o -96Hi,o,i - 128Hi,i,o - 96Hi jj) +4(3 -x)^3 

- S{l+2x-^ +x + 2x2)H_i,o + 8(l -x)(HiC2-Hi,o,o) -2(1 +x)(4H_i^2 + 4H0C2 

— 3Hi,o + 8H_i _i,o — 4H_i,o,o — 5Ho,o,o — 4Ho,o,i — 2Hoj,o — 2Hoj,i) +2(1 + 5x)Hi,i 

-^(3 + 79;c)Ho-8(4 + 5x + 2x2)^2 + ^(39-17x)Hi + (45 + 49;c+16;c2)Ho,o 

+ 8(3 + 5x)Hoj - ^ (233 - 131x) + 16Ho,-i,ox + { ^ " 78^3 + 69^2 + 6C2^}6( 1 -^)) 

+ C^Q (-8^3 + 2pqq{-x) (7^3 - 8H_i ^2 - 2Ho + 2H0C2 - 8H_i,_i,o + 10H_i,o,o 

+ 4H_ 1,0,1 + 6Ho,_i,o - 3Ho,o,o - 2Ho,o,i) - -^Pmi^) (3155 - 216^3 - 1584^2 + 4302Ho 
- 432H0C2 + 2202H1 - I296H1C2 + 1980Ho,o + 1584Ho,i + 1296Ho,_i,o + 648Ho,o,o 
+ 792Hi,o + 792Hi,i + 432Ho,o,i + 864Hi,o,o + 432Hi,o,i - 432Hi,i,o) - ^(3 - 5x)Hi 
+ 4(1+2x-'+x + 2x2)H_i,o-4(1-jc)(Hi^2-Hi,o,o)-4(4 + 3x + 2jc2)Ho,o 
+ 4(1 +x) (H_iC2 - 2Ho,i + 2H_i,_i,o - H_i,o,o) + 4(3 + 2x + 2x^%2 -\{29- 271x)Ho 

5 /ni ^ OTT r 5465 140^ 25K 71. 2U,, n\ 

+ — (91 + 167;c) - 8Ho,-i,ox+| " + — ^3 - —^2 + y ^2'|6(1 -x)) 

+ CpHf ^pqq {x) (247 - 1 44C2 + 342Ho + 1 74H 1 + 1 80Ho,o + 1 44Ho, 1 + 72H 1 ,0 + 72Hi , 1 ) 

+ ^(l-lU)Ho + (l-3x)Hi + ^(5-119x)+{^ + yC2 + ^C3}5(l-.^)) . (A.3) 
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Finally full three-loop coefficient function for p^'^^^ approximated by Eq. (13.71) is given by 

cW(^) = l_^i6Q(l-5x)HoAO,-i,o + y(l-2x)Ho,o,-iAO-^(l-^)(4376 

+ 3865Hi+684Ho,-i^3 + 1728Hi,o + 2730Hij-720Ho,-i,-i^2 + 432Ho,-i,o^2 

+ 318Hi,o,o - 198Hi,o,i + 198Hi,i,o + 216Hi,o,-i,o - 54Hi,o,i,o - 720Ho,-i,-i,-i,o 

+ 1044Ho,-i,-i,o,o + 360Ho,-i,-i,o,i + 504Ho,-i,o,-i,o - 504Ho,- 1,0,0,0 - 252Ho,-i,o,o,i 

2 

+ 72Ho,-i,o,i,o + 936Ho,o,-i,-i,o-72Ho,o,i,o,i+72Ho,o,i,i,o) + ^(l+^)(168Ho,iC3 

+ 12Ho,i,o^2 — 205Ho,i,i +60Ho,i,i^2 + 18Ho,i,o,i — 18Ho,i,i,o — 54Ho,o,o,i,o — 90Ho,o,o,i,i 

~6Ho,o,i,o,o — 24Ho,i,o,-i,o + 24Ho,i,o,o,o + 18Ho,i,o,o,i +6Ho,i,o,i,o — 54Ho,i, 1,0,0) 

4 2 10 

+ -(1 + 35x)Ho,o,-i ^2 - 3 (3 - 56x-^ + 3x - 56x2)H_i,_i,o,o - y (3 - 27x - 16x2)Ho,o,i,i 

40, , 80, 1 9, 

+ y(3-4x-i+3x-4x2)H_i,_i,_i,o-y(3 + 2x^i+3x + 2x2)H_i,o,i,i 

+ ^(3 + 4x-i-3x-4x2)(14Hi^3+5Hi,i^2)-y(3 + 7x-i+3x + 7x2)H_i,_i^2 

28 2 

- y (3 + 8x-i + 3x + 8x2)H_i,o,o,o - ^ (3 - 29x - 16x2)Ho,o,i,o + 8(4 + 3x)Ho,o,o^2 

- ^(6 - 7x-' - 6x + 7x2)Hi,o^2 - ^(6 - m - 35x2)Ho,o,o,i - ^(6 + 261x - 268x2)Ho,o,o 

- ^(7 + llx)Ho,o,-i,o,i - ^(9- 8x-i +9x- 8x2 )H_ 1,0,-1,0 _ ^(n +9x)Ho,o,o,o,i 

- y (9 + 8x-i + 9x + 8x2)H_i,o,o,i + ^ (9 - 8x - 26x2)Ho,oC2 + ^(11 + 15x)Ho,o,i ^2 
+ ^(15-4x-'-15x + 4x2)Hi,o,o,o-^(15 + 8x-i + 15x + 8x2)H_i,o,i,o 

+ ^(15 + 16x-i + 15x+ 16x2)H_i,o^2 + ^(15 - 3x- 8x2)Ho,o,-i,o " ^(28 + 15;c)C5 

- ^(21 +48X-1 - 9x + 68x2 )Ho,- 1,0,0 + ^(27 + 16x-i +9x- 8x2)Ho,iC2 _ 4^29 + 5x)^2C3 

8 4 

- — (27 + 134x"^ + 27x+ 134x2 )H_ 1,0,0 - -(33 - 16x"' - 15x- 8x2)Ho,_i,_i,o 

+ ^(33-4x-i-33x + 4x2)Hi,o,o,i + ^(33 + 39x + 28x2)Ho,i,o,o-^(41 + 36x)Ho,i,o 

4 2 

--(39+16x-^ -5x+16x2)Ho,_i,o,i --(39 + 28x-^ -39x-28x2)Hi,i,o,o 

2 4 

+ - (45 + 48x" ^ + 5x + 40x2 )Ho,_ 1 ^ 2 + ^ (5 1 + 64x" ^ + 5 Ix + 64x2 )H_ 1 1 ,0, 1 

+ ^(51 -9x- 110x2 )Ho^3 - ^(51 +20x)Ho^2' + ^(99- 150x-670x2)^3 

8 2 

- 27(m - 134x"^ + lllx- 134x2)H_i,_i,o + -(lll + 188x"' + lllx+ 188x2)H_i^3 

4 2 
-2^(111 +69x-268x2)Ho,o,i + 2^(123 + 268x^1 - 123x - 268x2)Hi^2 

+ ^(177 + 134x^1 + 177x+ 134x2)H_i^2 - ^^(186 -69x + 536x2)Ho^2 
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-^(221 - 18;c"' +221;c- 18;c2)H_i,o - ^(291 +201;c - 364;c2)^2^ 

4 8 

- — (297 -639.t + 268a;2)Ho,_i,o- — (321 + 134;*;-' + 321;c+ 134;c2)H_i,o,i 

2 2 2 

+ — (1172 - 1019;c+ 18;c2)^2 - — (1342 -77;c+ 18;c2)Ho,o - ^(1393 - 1019;c)Hoj 

- ^(4192 - 368U)Ho + ^(19Ho,oC3 - 6Ho,o,o,o,o) + yHo,o,o,ox' 

+ C| (-^;^qq(-^) (4960^5 + 150^3 - 140^2-7320^2^3 + 8I6C2' -4110H_i^3 

- 280H_i ^2 + 952H_ 1 ^2^ + 1225Ho - IOHq^s + 555Ho^2 + 772Ho^2^ + 30600H_i,_ 1 ^3 
+ 3960H_ 1,-1^2 + 240H_ 1 ,0 - 1 8600H_ 1 ,0 ^3 - 3 1 80H_ 1 ,0 - 22000Ho,_ 1 ^3 

- 2220Ho,_i ^2 + 550Ho,o + 5720Ho,oC3 - 1070Ho,o^2 + 260Ho,i - 3080Ho,i ^3 

- 1830Ho,i^2 - 30240H_i _i,_i^2 + 1920H_i,_i,o + 35760H_i _i,o^2 + 22560H_i,o,-i^2 

- 1580H_ 1,0,0 - 17880H_i,o,o^2 + 1240H_ 1,0,1 + 1640H_i,o,i^2 + 22680Ho,-i,-i^2 

- 1700Ho,-i,o - 23680Ho,-i,o^2 - 14000Ho,o,-i ^2 + 1 195Ho,o,o + 4520Ho,o,o^2 

- 1220Ho,o,i -3280Ho,o,i^2-240Ho,i,o-2120Ho,i,oC2-240Ho,i,i - 1720Ho,i,i^2 
+ 2160H_i,_i,_i,o - 1800H_i,_i,o,o - 2880H_i,_i,o,i - 900H_i,o,-i,o - 270H_i,o,o,o 
+ 3000H_i,o,o,i + 1440H_i,o,i,o + 1920H_i,o,i,i - 840Ho,-i,-i,o + 120Ho,-i,o,o 

+ 1800Ho,-i,o,i - 280Ho,o,-i,o + 1490Ho,o,o,o + 390Ho,o,o,i - 630Ho,o,i,o - 960Ho,o,i,i 
-60Ho,i,o,i+60Ho,i,i,o-5760H_i,_i,_i,_i,o + 25680H_i,_i,_i,o,o + 27360H_i,_,,_i,o,i 
+ 5520H_i,_i,o,-i,o-21800H_i,_i,o,o,o-31120H_i,_i,o,o,i -9920H_i,_i,o,i,o 

- 12480H_i,_i,o,i,i +5760H_ 1,0,- 1,-1,0 - 18800H_i,o,-i,o,o - 19680H_i,o,-i,o,i 

- 6560H 

1.0,0,-1,0 + 9480H_i,o,o,o,o + 15720H_i,o,o,o,i + 7000H_i, 0,0, 1, + 9040H_ 1,0,0, 1,1 
+ 3040H_i,o,i,o,o + 2160H_i,o,i,o,i + 1360H_i,o,i,i,o + 1440H_i,o,i,i,i+6000Ho,-i,-i,-i,o 

- 18840Ho,-i,-i,o,o-19680Ho,-i,-i,o,i-4800Ho,-i,o,-i,o + 13920Ho,-i,o,o,o 

+ 20120Ho,-i,o,o,i + 5760Ho,-i,o,i,o + 7280Ho,-i,o,i,i - 5920Ho,o,-i,-i,o + 1 1200Ho,o,-i,o,o 
+ 11040Ho,o,-i,o,i +3840Ho,o,o,-i,o -2040Ho,o,o,o,o -4200Ho,o,o,o,i -2440Ho,o,o,i,o 

- 3240Ho,o,o,i,i — 1280Ho,o,i,o,o — 1040Ho,o, 1,0,1 —720Ho,o, 1,1,0 —720Ho,o, 1,1,1 

- 160Ho,i,o,-!,o+480Ho,i,o,o,o + V20Ho,i,o,o,i - 160Ho,i,o,i,o + 480Ho,i,i,o,o + 80Ho,i,i,o,i 

-80Ho,i,i,i,o)-^Pqq(^)(15015 + 214400C5 -34320^3 + 35480^2 + 74240^2^3 

+ 15200Hi^2 + 48140Ho - 167040Ho^3 + 6488OH0C2 + 59520Ho^2^ - 311040Ho,-i,-i^2 

- 3136^2^ + 54528H1 ^2^ + 263680Ho,_ 1^3 + 130240Ho,_i^2 - 13480Ho,o - 35840Ho,o^3 

- 113280Ho,oC2 - 31640Ho,i + 148480Ho,i^3 -45120Ho,i^2 -47640Hi,o + 175360H 1,0^3 

- 84480Hi,o^2 - 33480Hi,i + 253440Hi,i ^3 - 51840Hi,i ^2 + 1 1220Hi - 17280Hi ^3 
+ 11840Ho,_i,o + 231680Ho,-i,o^2 + 220160Ho,o,-i^2 - 16240Ho,o,o - 35840Ho,o,o^2 
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- 86160Ho,o,i - 29440Ho,o,i^2 - 88800Ho,i,o - 14080Hoa,o^2 - 60160Ho,ia 

- 51200Ho,i,i ^2 - 55040Hi,o -1 ^2 - 44400Hi,o,o + 57600Hi,o,o^2 - 45920Hlo,i 

- 46080Hi,o,i ^2 - 59520Hi,i,o + 2560Hi,i,o^2 - 25920Hij,i - 46O8OH14 j ^2 

+ 120960Ho,-i,-Lo - 64640Ho,_ 1,0,0 - 69760Ho,- 1,0,1 - 34560Ho,o,-i,o + 30400Ho,o,o,o 
+ 91520Ho,o,o,i +28800Ho,o,i,o + 61440Ho,o,i,i - 6720Ho,i,o,o + 58560Ho,i,o,i 
+ 72000Ho,i,i,o + 77760Ho,i,i,i - 17280Hi,o,-i,o + 17280Hi,o,o,o + 90240Hi,o,o,i 
+ 76800Hi,o,i,o + 77760Hi,o,i,i + 2880Hi,i,o,o + 86400Hi,i,o,i + 103680Hi,i,i,o 
+ 86400Hi,i,i,i - 104960Ho,_i,_i,_i,o + 268800Ho,_i,_i,o,o+258560Ho,_i,_i,o,i 
+ 84480Ho,-i,o,-i,o - 130560Ho,- 1,0,0,0 - 186880Ho,-i,o,o,i - 23040Ho,-i,o,i,o 
-23040Ho,-i,o,i,i +276480Ho,o,-i,-i,o -230400Ho,o,-i,o,o -81920Ho,o,-i,o,i 

- 102400Ho,o,o,-i,o - 7680Ho,o,o,o,o +46080Ho,o,oAi + 84480Ho,o,o,i,o + 116480Ho,o,o,i,i 
+ 44800Ho,o,i,o,o + 1 24 160Ho,o, 1,0,1 + 142080Ho,o,i,i,o + 142080Ho,o, 1,1,1 

+ 74240Ho,i, 0,-1,0 - 5376OH0, 1,0,0,0 + 128000Ho,i,o,o,i + 139520Ho,i,o,i,o 

+ 140800Ho,i,o,i,i +75520Ho,i,i,o,o + 147200Ho,i,i,o,i + 172800Ho,i,i,i,o 

+ 138240Ho,i, 1,1,1 +243200Hi,o,-i,-i,o -48640Hi,o,-i,o,o + 176640Hi,o,_i,o,i 

-71680Hi,o,o,-i,o -65280Hi,o,o,o,o + 11520Hi, 0,0,0,1 + 131840Hi,o,o,i,o 

+ 135680Hi,o,o,i,i + 67840Hi,o,i,o,o + 138240Hi,o,i,o,i + 158720Hi,o,i,i,o 

+ 126720Hi,o,i,i,i + 120320Hi,i,o,_i,o -94720Hi,i,o,o,o + 145920Hi,i,o,o,i 

+ 161280Hi,i,o,i,o + 126720Hi,i,o,i,i + 88320Hi, 1,1,0,0 + 138240Hi,i, 1,0,1 

128 

+ 161280Hi,i,i,i,o + 1 15200Hi,i,i,i,i ) + 9^)Ho,o,o,-i,o - 72(1 - 3^)Hi,i,i,i 

+ 11^(1 _2^)Ho,o,-i,o,o -8(1 -^)(15Hi^2'-38Hi,o^3-6Hi,i^3 + 28Hi,o,-i^2 

- 60Hi,o,o^2 — 12Hi,o,i^2 — 56Hi, 1,0^2 — 12Hi, 1,1^2 — 8Hi,o,-i,-i,o — 12Hi,o,-i,o,o 
-32H 

1,0,-1,0,1 — 4Hi,o,o,-i,o + 20Hi,o,o,o,o + 48Hi,o,o,o,i +4Hi,o,o,i,o + 4Hi,o,o,i,i 

- 12Hi,o,i,o,o — I6H1, 1,0,-1,0 + 24Hi, 1,0,0,0 + 36Hi, 1,0,0,1 — 28Hi, 1,1,0,0) 

+ ^(1 +^)(177H_i^2' + 1530H_i,_i^3 - 630H_i,o^3 - 1860H_i,_i,_i^2 

+ 1440H_i,_i,o^2+1140H_i,o,-iC2-840H_i,o,o^2 + 60H_i,o,i^2-360H_i,_i,_i,_i,o 

+ 1500H_i,_i,_i,o,o + 1680H_i,_i,_i,o,i + 360H_i,_i,o,-i,o - 720H_i,_ 1,0,0,0 

- 1140H-i,-i,o,o,i - 120H- 1,-1, 0,1,0 - 120H-i,-i,o,i,i +360H-i,o,-i,-i,o 

- 900H_i,o,-i,o,o - 960H_ 1,0,- 1,0,1 - 300H_ 1,0,0,- 1,0 + 300H_i,o,o,o,o + 660H_i,o,o,o,i 
+ 60H_i 0,0,1,0 + 60H_ 1,0,0, 1,1 — 120H_ 1,0,1, 0,0 + 360Ho,o, 1,1,1 +250Ho,i,o,i,o 

+ 180Ho,i,o,i,i +210Ho,i,i,o,i +240Ho,i,i,i,o + 180Ho,i,i,i,i) - y (2 - 9^)Hi,o,i,o 

+ 16(3-1 l^)Ho,i,i ^2 - 64(4 + x)Ho,-i,o,-i,o - 8(5 - 39x)Ho, 1,0,0,1 - 8(5 - 27x)Ho, 1,0,0,0 
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16 32 
-y(5-ll;c)Ho,o,-iC2 + y(5 + ;c)Ho,-i,o,i,o + 64(5 + 2x)Ho,-i,-i,-i,o 

88 32 

- 4(5 + 57^)Ho,o,o^2 + 24(7 - 17x)Hoj,o^2 " y (7 - 13x)^2C3 " y (7 + 5x)Ho,o,-i,o,i 

32 64 64 

- y (9+ 16x)Ho,o,1^2 + y (10+ llx)Ho,o,l,l,0 - y (10 + 30x + 9x2)jj^^^ _^ ^ 

+ y ( 1 1 - x)Ho, 1 ,1 ,0,0 + y ( 14 + 7x- 1 - 44x - 7x2)Hi ,0^2 - y ( 17 + 79x)Ho,- 1 1 ,o, i 

8 16 32 

--(19-22x-i-75x + 22x2)Hi,o,o,i-y(19-47x)Hi,i,o,i-y(19-x)Ho,o,-i,-i,o 

+ 6(19 + 3x)Hi,i,i + 12(19 + 55x)Ho,i,i,i - ^ (23 - 97x)Hi,i,i,o - 4(23 - 45x)Ho,o,o,o,o 

+ - (25 - 24x"^ - X + 24x^)Hi,o,o,o - 32(25 + 2x"' + 25x + 2x2)H_i,o,i,i 

+ 16(27 +x)Ho,-i,o,o,o + y (29 + 41x)Ho,o,o,i,o + ^(31 - 361x)Ho,o^3 

-^(31-149x)Hi,o,i,i+y(31 + 33x)Ho,o,i,o,o + y(32 + 85x)Ho,-i,-i^2 

+ y (37 + 33x)Ho,o,i,o,i + y (43 + 83x)Ho,_i,o,o,i + y (44 + 1 16x + 1 lx2)Ho,i,o,o 

- ^(50 - 24x-i +259X+ 12x2)Ho,i^2 + ^(47 + 171x)Ho,o,o,o,i + ^(61 + 85x)Ho,o,o,i,i 

ao 1 ^ 

- 8(43 + 87x)Ho,-iC3 " y (^5 + 6x-^ +55x + 6x2)H_i,o,i,o - y (67 + 101x)Ho,-i,o^2 

32 16 

+ y (64 + 78x"' +73x + 78x2)H_i,_i,o,i + y (§1 - 12^"' +43x- 12x^)H_i,_i,_i,o 

_ ^(95 - 48x-i - 33x + 48x2)Hi,o,-i,o + y (68 + 59x)Ho,i,o - y (95 + 91x)Ho,-i,-i,o,o 

8 8 
--(113-24X"' +59x-24x2)H_i,o,-i,o + ^(115 + 216x"^ +203x + 216x2)H_i,_i,o,o 

+ -(119+ 12x-' - 137x-12x2)Hi,i^2 + ^(122-65x)Hi,o,i -y(103 + 120x)^5 

-^(145-48x"' -95x-24x2)Ho,_i,_i,o + ^(153 + 116x"^ -281x- 116x2)Hi^3 

- ^ ( 165 + 48x- ' + 253x + 2 1 6x2 )Ho,- 1 ^^^^ _ ^ ^ ^75 ^ 324;,-- 1 + 249x + 324x2 )H_ i,_ 1 ^2 
+ A(i29 - 151x)HoC2' + ^(199 + 495x)Ho,i,i,o - -^(203 - 720x - 148x^)^2^ 

- ^ (227 + 44x- ' - 269x - 44x2 ^ jj^ ^ - ^ (26 1 - 48 8x- ^ + 23 Ix - 48 8x2 )H_ 1 _ 1 ,0 
+ ^ (204 + 463x + 24x2 )Ho,o,i,o + ^ (253 + 597x + 24x2)Ho,o,i,i + ^ (101 + 29x)Hi,i,o 

- -(267 + 488x^1 +276x + 488x2)H_i,o,o - y (275 + 108x^1 +284x+ 108x2)H_i,o,o,i 

8 20 4 

+ -(294 + 264X-1 + 359x + 264x2 )H_ 1^3 + y (394 - 427x)Hi + -(245 + 61 3x)Ho, 1,0,1 

-^(455 + 48x^1 + lllx + 312x2)Ho,-i,o,i -^(471 -488x"i - 159x + 488x2)Hi^2 

+ ^ (53 1 + 488x^1 + 541x + 488x2)H_ 1 ^2 + ^ (541 + 252x-i + 555x + 252x2)H_ 1,0^2 
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- ^(599+ 192x^1 + 641a; + 192x2 )H_ 1,0,0,0 - ^(669- 1629a; + 97 6a;2)Ho,_ 1,0 

+ ^(755 + 1383.x + 480a;2)Ho,o,o,o + ^(765 + 144a;-' + 317a; + 648a;2)Ho,- 1^2 

12 8 

+ -(817-763a;)Hi,o,o + -(851+615a;)Ho,i,i --(927 + 488A;-' +927a; + 488x2)H_i,o,i 

4 2 

- - ( 1045 + 462.x- 1 + 877a; + 1 S6x^)n-i ,0 + - ( 1 975 + 2847a; + 776a;2 )Ho,o,o, 1 

-^(2035 + 2271a;+ 1288a;2)Ho^3 - ^ (2295 + 2847a; + 1064a;2)Ho,o^2 

+ I- (2823 + 1305a; + 3904a;2 )Ho + 7 (4143 - 5305a;)Hi - ^ (6705 + 3537a; + 9760a;2 )^3 
9 '6 '9 

+ ^ (7533 + 6969x + 3904a;2 )Ho,o,i - ^ ( 10209 + 6969a; + 7808a;2 )Ho^2 - 128Ho,i,o,- 1,0^ 

+ ^ (7129 - 741 lx)Hi 1 + -!_ ( 1 1 17 1 _ 29353a; + 1488a;2 )Ho + A ( 12 137 - 30523;c)Ho 1 
6 18 18 

- ^(20497 - 30523a; + 1488;c2)^2 + ^(28579 + 7263a;)Ho + -^(481783 -463765a;) 

4 4 f 7255 950 

+ 16(17Ho,i^3 +4Ho,-i,o,u) + 3^1 (x) + -g2{x) + | + 1240^5 - —^3 

304 2 1129 2148 2^4184 3^ 

-— ^3 2~^2 + 808Q2Q3 ^Q2 +^75-^2 >0{l-x] 

+ Cjnf (^Apqq(_;c) (2700^3 + 830^2 + 1026^2^ - 1440H_i^3 - 3000H_i^2 - 1305Ho 
+ 990Ho^3 + 345H0C2 + 1 800H_ 1,-1^2 + 1 660H_ 1 ,0 - 1260H_ 1 ,0^2 - 2340Ho,_ 1 ^2 

- 1910Ho,o + 450Ho,o^2 - 3600H_i,_i,o + 4800H_i,o,o + 1200H_i,o,i + 3000Ho,_i,o 

- 1245Ho,o,o - 600Ho,o,i + 3600H_i,_i,_i,o - 6120H_i,_i,o,o - 3960H_i,o,-i,o 
+ 6 1 20H_ 1 ,0,0,0 + 360H_ 1 ,0,0, 1 - 720H_ 1 ,0, 1 , 1 - 3960Ho,_ 1 1 ,0 + 6 1 20Ho,_ 1 ,0,0 
+ 360Ho,- 1,0,1 +4500Ho,o,-i,o -2070Ho,o,o,o -360Ho,o,o,i +360Ho,o,i,i) 

^ -;7qq(x) (30045 + 580080^3 +92280^2 - 89568^2^ + 287700Ho + 420480Ho^3 
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+ 505680Ho^2 + I494OH1 + 187200Hi^3 + 253440Hi^2 - 138240Ho,_i^2 + 64920Ho,o 
+ 262080Ho,o^2 - 145400Ho,i +92160Ho,i^2 - 306680Hi,o + 43200Hi,o^2 -245880Hi,i 
+ 57600Hi,i^2 + 224640Ho,_i,o - 196080Ho,o,o - 489360Ho,o,i - 333120Ho,i,o 

- 371520Ho,i,i - 301200Hi,o,o - 368640Hi,o,i - 387360Hi,i,o - 302400Hi,i,i 

+ 276480Ho,_ 1,0,0 + 138240Ho,-i,o,i +276480Ho,o,-i,o - 129600Ho,o,o,o -267840Ho,o,o,i 

- 279360Ho,o, 1,0 -2908 80Ho,o, 1,1 -204480Ho,i,o,o -218880Ho,i,o,i -293760Ho, 1,1,0 

- 236160Ho,i 1.1 + 138240Hi.o -1,0 - 152640Hi,o,o.o - 100800Hi,o.o 1 - 264960Hi 0,1,0 



- 195840Hi,o,i,i - 262080Hi,i,o,o - 172800Hi, 1,0,1 - 236160Hi,i,i,o - 172800Hi,i,i,i ) 
+ ^(1 -39x+ 14x2)Ho,i,o + ^(1 -7x)Ho,i,o,o - y (1 -x)(3Hi^3 + 7Hi,o^2 + 5Hi,i^2 

— 4Ho,-i,o,i — 4Hi,o,-i,o — 6Hi,o,o,o — 5Hi,o,o,i + Hi,o,i,o — Hi, 1,0,1 + Hi, 1,1,0) 
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+ 1 ( 1 + ^) (48H_ 1 ^3 - 48H_ 1 1 ^2 + 24H_ 1 ,0 ^ 2 + 69Ho,o ^2 - 96H_ i i i ,o 

+ 1 44H_ 1 1 ,0,0 + 96H_ 1 ,0,- 1 ,0 - V2H_ 1 ,0,0,0 + 96Ho,- 1 1 ,0 - 9 1 Ho,o,o,o - 69Ho,o,o, 1 

-42Ho,o,i,o -42Ho,o,i,i - I8H0, 1,0,1 - I8H0, 1,1,0 - 18Ho,i,i,i) - y (1 +7^)Ho,- 1,0,0 

- ^{l+2lx+ 14x2)Hi,o,i + ^{2-x-^ +2x-x^)}i^i^^ifl + ^{591 - 1133;c)Hi,i 

- ^(3 - Ux)Ho^3 + ^(5 - 117x)Ho,i,i -^(7 + 4x-' -29x- 18^2)jjj^2 

64 64 

+ —{S + 3x-^+Sx + 3x^)U^ifl^i + —{S + 9x-^ + Sx + 9x^)U^iflfl 

32 1 6 

-— (10 + 9;c"' + 10;c + 9A;2)H_i^2 + y(15 + 24;c"i - llx+ 12;c2)Ho,_i,o 

8 4 16 2 

- -(5 - llx)Ho,i^2 - — (23 + 63x)^2' - y (37 + 43x+ 12x2)Ho,o,o + 3 (59 - 113x)Hi,o,o 

+ ^(75 - 106x)Hi + ^(103 - 133x + 84x2)Hi,i,o - ^(108 + 221x + 66x^)Ho,o,i 
+ ^(138 + 221x + 90x2)Ho^2 + ^(158 + 141x"' + 158x+141x2)H_i,o 

- ^(129 - 325X+ 132x2)^3 - i^(341 + 5617x)Ho,i + y (2Ho,-i^2 - 3Ho,o,-i,o)^ 

+ (625 - 3283x)Hi,o + ^ (2869 + 5617x + 2256x2 )^2 + ^ (4105 - 291 lx)Ho 
27 ' 27 81 

1 1 32 

- —(16627 + 13429X + 6768x2 )Ho + —(32905 - 2097 Ix) + — (2 - 5x)Hi 1 1 
81 216 3 

r 341 592^ 1348^ 549K 352^^ 16472^ 2U/, ^ 
+ - + - y^^2 - ^^2^3 - j5(l -X, 

+ Cpnf (^-^p^^{x){4351 - 864^3 - 7236^2 + II6IOH0 - 3888Ho^2 + 4230Hi 

- 1296Hi^2 + 10980Ho,o + 7236Ho,i +3132Hi,o + 3132Hi,i +4968Ho,o,o + 3888Ho,o,i 
+ 2592Ho,i,o + 2592Ho,i,i + 1296Hi,o,o + 1296Hi,o,i + 1296Hi,i,o + 1296Hi,i,i ) 

+ ^(1 - 3x)(10C2 - 19Hi - 10Ho,i - 6Hi,o - 6Hi,i) - ±(17 -59x)Ho,o - ^(171 -577x) 

2 r 9517 152. 860. 32. 2 



- -(197 - 803.)H„ _C/|8(1 

+ C|Q (^4^7?qq(-^)(197640^5 - 98550^3 - 92260^2 - 338040^2^3 + 33102^2^ 

- 47610H_i ^3 + 97080H_i ^2 + 35856H_i ^2^ + 139905Ho - 99990Ho^3 - 17655Ho^2 
+ 38124Ho^2^ + 1608120H_i,_i^3 + 27720H_i,_i^2 - 78680H_i,o - 860760H_i,oC3 

- 33660H_i,o^2 - 1088640Ho,_i^3 + 43020Ho,_i^2 + 137620Ho,o + 232200Ho,o^3 

- 133470Ho,o^2 + 52920Ho,i - 182520Ho,i^3 - 107730Ho,i^2 - 1594080H_i,_i,_i^2 
+ 244800H_ 1 1,0 + 1 864080H_ i,_ 1 ,0 ^2 + 1 1 34000H_ 1 ,0,- 1 C2 - 276 1 80H_ 1,0,0 

- 854280H_ 1 ,0,0 C2 + 25320H_ 1 ,0,1 + 1 1 3400H_ 1 ,0,1 ^2 + 1 143720Ho,_ 1,-1^2 

- 208860Ho,_i,o - 1 185840Ho,_i,o^2 - 654480Ho,o,-i C2 + 143385Ho,o,o + 191 160Ho,o,o^2 
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-28860Ho,o,i - 183600Ho,o,i^2-6480Ho,i,o- 126360Ho,i,o^2-6480Ho,i,i 

- 98280Ho,i,i^2 - 100080H_i _i _i,o + 220680H_i_ 1,0,0 - 77760H_i,_ 1,0,1 

+ 149940H_i, 0,-1,0 - 273330H_i,o,o,o + 65160H_i,o,o,i + 38880H_i,o,i,o + 83520H_i,o,i,i 

+ 151560Ho,-i,-i,o-272520Ho,-i,o,o + 32760Ho,-i,o,i- 183960Ho,o,-i,o 

+ 173970Ho,oAO + 84690Ho,oAi -17010Ho,o,i,o-41760Ho,o,i,i+6480Ho,i,o,o 

- 1620Ho.i,o.i + 1620Ho,i,i,o - 155520H_i,_i,_i,_i,o + 1 134000H_i,_i,_i,o,o 

+ 1516320H_i,_i,_i,o,i + 123120H_i,_i,o,-i,o - 916920H_i,_i,o,o,o - 1643760H_i,_i,o,o,i 

- 509760H_i,_i,o,i,o - 682560H_i,_i,o,i,i + 142560H_i,o,-i,-i,o - 792720H_ 1,0,- 1,0,0 

- 1062720H_i,o,-i,o,i - 194400H_i,o,o,-i,o + 381240H_i,o,o,o,o + 735480H_i,o,o,o,i 
+ 294840H_i,o,o,i,o +416880H_i,o,o,i,i + 1425 60H_ 1,0,1, 0,0 +79920H_ 1,0, 1,0,1 

+ 1 5 120H_ 1 ,0, 1 , 1 ,0 + 38880H_ 1 ,0, 1 , 1 , 1 + 1 62OOOH0,- 1 1 1 ,0 - 79596OH0,- 1 1 ,0,0 

- 1062720Ho,-i,-i,o,i - 120960Ho,-i,o,-i,o +568080Ho,-i,o,o,o + 1011960Ho,-i,o,o,i 
+ 280800Ho,-i,o,i,o + 369360Ho,-i,o,i,i - 142560Ho,o,-i,-i,o +453600Ho,o,- 1,0,0 

+ 583200Ho,o,-i,o,i + 125280Ho,o,o,-i,o - 81000Ho,o,o,o,o - 178200Ho,o,o,o,i 

- 91800Ho,o,o,i,o - 130680Ho,o,o,i,i - 56160Ho,o,i,o,o - 36720Ho,o, 1,0,1 - 10800Ho,o, 1,1,0 

- 19440Ho,o,i,i,i - 12960Ho,i,o,-i,o + 30240Ho,i,o,o,o +49680Ho,i,o,o,i - 12960Ho,i,o,i,o 
+ 30240H0, 1 , 1 ,0,0 + 6480Ho, 1 , 1 ,0, 1 - 6480Ho, 1 ,1 , 1 ,0 ) + ^7?qq {x) (2292435 + 48 12480^5 
+ 2680080^3 + 1050360^2 - 1347840^2^3 - 1746432^2 V 5 132100Ho - 1625760Ho^3 
+ 4330560Ho^2 + 27 1296Ho^2^ + 1001340Hi - 4821 120Hi,o,-i ^2 + 25920Ho,_i,o,o 

+ 1565568Hi^2V6048000Ho,_i^3 +2967840Ho,_iC2 + 1045320Ho,o + 2056320Ho,o^3 

- 188640Ho,o^2 - 1679800Ho,i +7568640Ho,iC3 + 1260000Ho,i^2 - 3855640Hi,o 

+ 9849600Hi,o^3 + 501120Hi,o^2 -3033000Hi,i + 10108800Hi,i^3 + 1229760Hi,i^2 

- 7724160Ho,_ 1,-1^2 + 1 853280Ho,_ 1,0 + 437 1 840Ho,_ 1,0 C2 + 2695680Ho,o,- 1 ^ 2 

- 851520Ho,o,o + 1192320Ho,o,o^2 - 5074800Ho,o,i + 1762560Ho,o,i ^2 -4427520Ho,i,o 
+ 1866240Ho,i,o^2 -4648320Ho,i,i + 1503360Ho, 1,1^2 + 1838880Hi^3 + 2321280Hi^2 

+ 4008960Hi,o,o^2 -4279680Hi,o,i + 1520640Hi,o,i^2 - 5107680Hi,i,o + 1831680Hi, 1,0^2 
- 3741120Hi,i,i + 1313280Hi,i,i^2 + 3499200Ho,_i,_i,o - 3673200Hi,o,o 

- 1218240Ho,_i,o,i + 1382400Ho,o,-i,o -99360Ho,o,o,o -591840Ho,o,o,i -2904480Ho,o,i,o 
-2888640Ho,o,i,i -2897280Ho,i,o,o -2472480Ho,i,o,i - 3114720Ho,i,i,o -2597760Ho,i,i,i 
+ 1391040Hi,o,-i,o -2158560Hi,o,o,o - 1005120Hi,o,o,i -2914560Hi,o,i,o 

- 3841920Hi, 1,0,0 -2030400Hi, 1,0,1 -2364480Hi, 1,1,0 - 1900800Hi,i,i,i -2257920Hi,o,i,i 
-2937600Ho,_i,_i,_i,o + 6462720Ho,_i,_i,o,o + 6255360Ho,_i,_i,o,i +2246400Ho,_i,o,-i,o 
-2453760Ho,-i,o,o,o -3144960Ho,-i,o,o,i +241920Ho,-i,o,i,o + 380160Ho,-i,o,i,i 
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+ 7050240Ho,o,-i,-i,o -3525120Ho,o,-i,o,o + 829440Ho,o,-i,o,i " 1589760Ho,o,o,-i,o 

- 1296000Ho,oAO,o -984960Ho,o,o,o,i -535680Ho,oAi,o -535680Ho,o,o,i,i 

- 1866240Ho,o,i,o,o - 622080Ho,o, 1,0,1 - 138240Ho,o,i,i,o - 311040Ho,o,i,i,i 

+ 3628800Ho,i,o,-i,o -4320000Ho,i,o,o,o +915840Ho,i,o,o,i - 120960Ho,i,o,i,o 

-328320Ho,i,o,i,i - 2488320Ho,i, 1,0,0 -207360Ho,i,i,o,i +535680Ho,i,i,i,o 

+ 6186240Hi,o,-i,-i,o + 1900800Hi,o,-i,o,o+7914240Hi,o,-i,o,i+691200Hi,o,o,-i,o 

-4337280Hi,o,o,o,o -2108160Hi,o,o,o,i +241920Hi,o,o,i,o +414720Hi,o,o,i,i 

-2661120Hi,o,i,o,o -380160Hi,o,i,o,i + 103680Hi,o,i,i,o -311040Hi,o,i,i,i 

+ 5O8O32OH1 J, 0,-1,0 - 6272640H 1,1, 0,0,0 + 1693440Hi, 1,0,0,1 + 103680Hi,i,o,i,o 

- 241920Hi,i,o.i,i - 2937600Hi j,i,o,o - 69120Hi,i j,o,i + 622080Hi,i, 1,1,0) 

- 4(1 + 8x"' - 26x - 4x'^)Uo,\ ^2 - ^ (1 + 40x"' -2x- 40x^)Hi,o,o,i - 4(1 - 3x)Hi, 0,1,0 
320 

- — ( 1 - 2x)Ho,o,- 1 ,0,0 + 4( 1 - x) (7Hi ^2' - 94Hi ,0^3 - 22Hi , 1 ^3 + 92Hi ,0,- 1 ^2 

- 108Hi,o,o^2 — 12Hi,o,i^2 — 72Hi, 1,0^2 — 12Hi, 1,1^2 — 8Hi,o,-i,-i,o — 44Hi,o,-i,o,o 

- 96Hi,o,-i,o,i — 20Hi,o,o,-i,o + 36Hi,o,o,o,o + 80Hi,o,o,o,i +4Hi,o,o,i,o + 4Hi,o,o,i,i 

- 12Hi,o,i,o,o — 48Hi, 1,0,-1,0 +4OH1, 1,0,0,0 + 36Hi, 1,0,0,1 — 28Hi, 1,1,0,0) 

- ^(1 +x)(153H_i^2' + 2970H_i,_i^3 - 1350H_i,oCs - 3780H_i,_i,_i^2 

+ 2640H_i,_i,o^2 + 2100H_i,o,-i^2- 1680H_i,o,o^2 + 60H_i,o,i^2 + 195Ho,i,i,i 
-360H_i,_i,_i,_i,o + 2460H_i,_i,_i,o,o + 3600H_i,_i,_i,o,i+360H_i,_i,o,-i,o 

- 960H_i,_ 1,0,0,0 - 2100H_i,_i,o,o,i - 120H_i,_i,o,i,o - 120H_i,_i,o,i,i + 360H_i,o,-i,-i,o 

- 1380H_ 1,0,- 1,0,0 - 1920H_i,o,-i,o,i - 540H_ 1,0,0,- 1,0 + 540H_ 1,0,0,0,0 + 1260H_i,o,o,o,i 

+ 60H_i,o,o,i,o + 60H_i,o,o,i,i-240H_i,o,i,o,o+40Ho,i,o,i,o-60Ho,i,o,i,i+60Ho,i,i,i,o) 

8 32 16 32 

+ -(1 +2x)Hi,i,o,i + y (1 + 19x)Ho,i,o,-i,o + y (2 - 23x)Ho,o,-i^2 - y (2 + x)Ho,-i,o,i,o 

+ ^(3 + llx)(Ho,o,i,o,i - Ho,o,i,i,o) - ^(5 +4x-i - 9x -4x2)Hi,oC2 + ^(19 - 35x)Hi,i,i,o 

16 32 8 

- y (5 -41x)Ho,o,o,-i,o + y ~ 13x)Ho, 1,0,0,0 + j(5 + 3-^)(6Ho,o,o,i,o +Ho,o,i,o,o) 

16 8 32 

+ y (7 -26x)Ho,i,o,o,i - ^(7 + 5x)(5Ho,o,i^2 -6Ho,o,o,i,i) + y (7 + 8x)Ho,o,-i,o,i 

4 32 64 

--(21-31x)Hi,o,i,i-24(9 + 4x)Ho,o,o^2 + y(ll + 85x)Ho,-i,-i,o,i-y(ll-6x)Ho,o^3 

-16(14 + 57x)Ho,-i,-i^2 + y(16-7x)Ho,o,-i,-i,o-^(121 + 362x + 24x2)Ho,o,i,o 

16 , , 16 , , 16 , ^ 

+ y (8 + 23x)Ho,i , 1 ,0,0 - y (25 + 86x)Ho,- 1 ,0,0, 1 + y (3 1 - x)Ho,- 1 ,0,- 1 ,0 

-4(37 + 4x-i-23x-4x2)Hi,i^2 + y(37-10x)^2C3-^(43-89x)Ho,i,o^2 

32 4 16 

+ Y^(43 + 45x)Ho^2^ + ^(59 + 80x-i - 13x- 80x2)Hi,i,o,o + y (62 + 79x)Ho,-i,-i,o,o 
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+ ^(71 + 36^)Ho,oAO,i + ^(79 + 209^)Ho,- 1,0^2 + ^(^21 - 24^"' + 19x - 24^2)h_i „ _i q 

- ^(80 - 137^)Hi,i,i - ^(95 + 13^)Ho,- 1,0,0,0 + ^(107 - 144^-^ - 53^+ 144;c2)Hi, 0,-1,0 

8 8 
+ -(107 + 147;c + 84;c2)Ho,o,-i,o - ^ (109 - llx'^ +73;c- 12;c2)H_i,_i,_i,o 

16 4 
+ y(127 + 6;c-' + 127;c + 6;c2)H_i,o,i,o--(140+ 164;c"i -95;c- \6Ax^)\iit,^ 



-4(165 + 136^;"^ + 187;c+ 136;c2)H_i,_i,o,o - ^(193 + 617;c + 24;c2)Ho,o,i,i 



16 16 
+ y (175 + 6;c"^ + 175a; + 6;c2)H_i,o,i,i - y (200+ 174;c"i +209;c+ 174;c^)H_i,_i,o,i 

4 T 4 o 32 

- -(197 + 15;c + 80x2 )Ho, 1,0,0 _ _(237 -408x- 154x2)Hi,o,i - y (20 + x)Ho,_i,_i,_i,o 

- ^ (2 1 1 + 3 1 lx)Ho, 1 ,0, 1 + ^ (379 + 94x- ^ + 406x + 94^^ )H_ 1 ,0,0 - ^ (227 + 77x)Ho, 1 ^3 

- ^(273 - 216x + 154x2)Ho,i,o + ^(217 - 48x"^ - 155x - 24x2)Ho,_i,_i,o 

- 4(429 + 164x^1 + 433x + 164x2)H_i,o^2 - (455 - 661x)Hi,o,o + 2(269 + 201x)^5 

+ ^ (465 - 622x^1 - 1581x - 302x2)Hi ^2 + ^ (549 - 1056x"' - 3917x + 1244x2)Ho,_i,o 

2 2 

- - (563 - 144x- ' - 58 Ix + 144x2 )Hi ,0,0,0 ^ _ ( ^743 _^ 5 ^48^^ _ ^ ^^^2 

+ ^ (647 + 204x^1 + 656x + 204x2 )H_ 1,0,0,1 + ^(689 + 48x"^ + 289x + 408x2 )Ho,_ 1,0,0 
+ ^ (69 1 + 708x- ^ + 763x + 708x2 )H_i,_ 1^2-^ (898 + 552X-' + 96 Ix + 552x2 )H_ 1^3 

- ^(929 - 1415X + 924x2)Hi,i,o + ^(957 - 622x"i + 867x - 622t2)H_i,_i,o 

+ ^ ( 1 1 87 + 288x" ' + 1 205x + 288x2 )H_ 1,0,0,0 - ^ ( 1735 + 223x + 1296x2 )Ho,o,o,o 
-^(1313 + 1469x + 656x2 q q 1 + 1 ( 135 1 + 48,x" ' + 75x + 696x2) Ho,- 1,0,1 

+ ^(1617 + 192x + 796x2)^2^-l(1270-633x)Ho,i,i-^(21-71x)Ho,i,i,o 

2 8 16 

+ - (1741 + 1469X + 992x2 )Ho,oC2 _ _(582 + 145x + 80x2)Ho,o,i - y (20 -x)Ho,i,i^2 

+ ^(2281 +559X+ 1432x2)Ho^3 - ^(2485 + 144x-' +305x+ 1416x2)Ho,-i^2 

4 2 
+ - (2503 + 622x" ' + 2503x + 622x2 )H_ 1 ,0, 1 + - (2877 + 5 80x + 1 564x2 ^^^^^ 

+ ^(3219 + 1714X + 5882x2)^3 - ^ (4049 + 1 866x- ^ +4139x+ 1 866x2 )H_ 1^2 

2 1 

-—(9253 + 3882x^1 +9973x + 4926x2)H_i,o + —(28121 +235007x + 591 12x2 )Ho,o 
27 '162 

- "1^(29169 - 34843x)Hi,i +^-(40879 - 108101x - 19704x2 )^2 + ^(167 + 51 lx)Ho,_i ^3 
18 ' 54 3 

- (49721 - 81983x)Hi,o - ^(50101 - 52955x)Hi - ^(77891 - 108101x)Ho,i 
54 ' 18 54 ' 
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- ^(80219 - 7801 U)Ho - -^-(1449755 - 122560U) - 8(4Ho,-i,o,i.i - 21Ho,o,o,o,o) 
ol 432 

4,4^^ [9161 4952^ ^^^^^ 1016^ 2 104189^ 
6644^ „ 87632 „ 2 33556 „ 3,^,, 

-^^2^3 + -^C2 -^^2 )5(1-.) 

+ CpC^ (-^/'qql-^) (3 1860^5 - 51300^3 - 44240^2 - 70200^2^3 + 5535^2^ 

+ 31680H_i^3 + 52320H_i^2 + 5076H_i^2^ + 53415Ho-49860Ho^3- 16320Ho^2 
+ 8640Ho^2^ + 390960H_ 1,-1^3- 39600H_ 1,-1^2- 42580H_ 1 ,0 - 179280H_ 1 ,0 ^3 
+ 26100H_i,o^2 - 247320Ho,_iC3 + 51480Ho -1^2 + 61385Ho,o + 38880Ho,o^3 

- 52290Ho,o^2 + 2295OH0J - 49680Ho,i ^3 - 29160Ho,i ^2 - 388800H_i,_i,_i ^2 
+ 96480H_ 1 _ 1 ,0 + 449280H_ 1 1,0 ^ 2 + 262440H_ 1 ,0,- 1 C2 - 1 1 67 60H_ 1,0,0 

- 1 85760H_ 1 ,0,0 C2 - 4080H_ 1 ,0, 1 + 34560H_ 1 ,0, 1 C2 + 265680Ho,_ 1,-1^2- 8 1480Ho _ 1,0 

- 273240Ho,-i,o^2 - 138240Ho,o,-i ^2 + 55560Ho,o,o + 34560Ho,o,o^2 + 2040Ho,o,i 
-47520Ho,o,i^2-34560Ho,i,o^2-25920Ho,ij^2-79200H_i,_i _i,o + 134640H_i,_i,o,o 
+ 87120H_i,o,-i,o - 133020H_L0A0 -7920H_i,oAi + 15840H_Lo,Li +87120Ho,-L-i,o 

- 137880Ho,-i,o,o - 7920Ho,-i,o,i - 88200Ho,o,-lo + 66870Ho,o,o,o + 37080Ho,o,o,! 

- 7920Ho,o,Li + 3240Ho,Lo,o + 220320H_i,_l_ 1,0,0 + 388800H_l_i,_lo,i 

- 12960H_i,_Lo_i,o - 164160H_i,_Lo,o,o -401760H_i,_i,o,o,i - 120960H_i _i,o,i,o 

- 172800H_i _i,o,i,i -6480H_i,o,-i,-i,o- 142560H_ 1,0,- 1,0,0 -265680H_ 1,0,-1, 0,1 

- 8640H_i,o,o,-i,o + 62640H_ 1,0,0,0,0 + 155520H_i,o,o,o,i + 52920H_ 1,0,0, 1,0 
+ 86400H_i,o,o,i,i +30240H_i,o,i,o,o + 10800H_i,o,i,o,i - 10800H_i,o,i,i,o 

- 143640Ho,-i,-i,o,o -265680Ho,-i,-i,o,i +4320Ho,-i,o,-i,o +96120Ho,-i,o,o,o 
+ 234360Ho,-i,o,o,i +62640Ho,-i,o,i,o + 86400Ho,-i,o,i,i +8640Ho,o,-i,-i,o 

+ 75600Ho,o,- 1,0,0 + 142560Ho,o,-i,o,i + 10800Ho,o,o,-i,o - 12960Ho,o,o,o,o 
-32400Ho,o,o,o,i - 12960Ho,o,o,i,o - 21600Ho,o,o,i,i - 1 0800Ho,o,i, 0,0 -4320Ho,o, 1,0,1 
+ 4320Ho,o,i,i,o-4320Ho,i,o,-i,o + 8640Ho,i,o,o,o + 15120Ho,i,o,o,i -4320Ho,i,o,i,o 

+ 8640Ho,i,i,o,o+2160Ho,i,i,o,i-2160Ho,i,i,i,o)-:^;?qqW(2996875 + 65 1240^5 
-797580^3 -3438810^2 -252720^2^3 -20412^2 V5966955Ho - 1296000Ho^3 

- 1808460Ho^2 + 252720Ho^2V2281005Hi -324000Hi^3 - 1694520Hi^2 

+ 1 183896Hi ^2^ + 1399680Ho,_i ^3 + 680400Ho,_i ^2 + 4599450Ho,o + 699840Ho,oC3 

- 1607040Ho,o^2 + 3055590Ho,i + 2099520Ho,i ^3 - 1 125900Ho,i ^2 + 1062450Hi,o 
+ 3285360Hi,oC3 - 1474200Hi,o^2 + 1255230Hi,i +2799360Hi,i^3 - 891000Hi,i^2 

- 1982880Ho,_i,_i ^2 + 952560Ho,_i,o + 699840Ho,_i,o^2 - 155520Ho,o,_i ^2 
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+ 2502900Ho,o,o -V7760Ho,o,o^2 + 1551420Ho,o,i -77760Ho,o,i^2 + 1070820Ho,i,o 

- 213840Ho,i,o^2 + 784080Hoj,i + 38880Ho,i,i ^2 - 2293920Hi,o -1 ^2 + 887760Hi,o,o 
+ 252720Hi, 0,0^2 + 826200Hi,o,i - 97200Hi,o,i ^2 - 42120Hi,i,o - 349920Hi, 1,0^2 

+ 392040Hi,ij +77760Hi,i,i^2 + 1049760Ho,_i,_i,o + 505440Ho -1,0,0 - 155520Ho,-i,o,i 
+ 1040040Ho,o,-i,o + 1203660Ho,o,o,o + 1333260Ho,o,o,i + 306180Ho,o,i,o + 362880Ho,o,i,i 
+ 220320Ho,i,o,o + 341820Ho,i,o,i - 312660Ho,i,i,o + 913680Hi,o,-i,o + 541080Hi,o,o,o 
+ 1 146960Hi,o,o,i - 291600Hi,o,i,o + 320760Hi,o,i,i - 463320Hi, 1,0,0 + 178200Hi, 1,0,1 
-498960Hi,i,i,o - 855360Ho,_i,_i,_i,o + 1594080Ho,_i,_i,o,o + 1555200Ho,_i,_i,o,i 
+ 622080Ho,-i,o,-i,o - 388800Ho,-i,o,o,o - 349920Ho,-i,o,o,i +311040Ho,-i,o,i,o 
+ 388800Ho,-i,o,i,i + 1866240Ho,o,-i,-i,o -233280Ho,o,-i,o,o + 1088640Ho,o,-i,o,i 

- 116640Ho,o,o,-i,o -233280Ho,o,o,o,o + 116640Ho,o,o,o,i + 116640Ho,o,o,i,o 
+ 194400Ho,o,o,i,i -427680Ho,o, 1,0,0 + 1477440Ho,i,o,-i,o - 1263600Ho,i,o,o,o 

+ 913680Ho,i,o,o,i - 58320Ho,i,o,i,o + 19440Ho,i,o,i,i - 1088640Ho,i, 1,0,0 - 38880Ho,i,i,o,i 

+ 19440Ho,i,i,i,o + 1632960Hi,o,-i,-i,o + 1438560Hi,o,-i,o,o + 3110400Hi,o,-i,o,i 

+ 933120Hi,o,o,-i,o - 1127520Hi,o,o,o,o +291600Hi,o,o,o,i +505440Hi,o,o,i,o 

+ 777600Hi,o,o,i,i - 1049760Hi,o,i,o,o + 38880Hi,o,i,o,i -38880Hi,o,i,i,o 

+ 1944OOOH1, 1,0,-1,0 -2021760Hi,i,o,o,o + 1205280Hi, 1,0,0,1 - 116640Hi,i,o,i,o 

+ 77760Hi,i,o,i,i - 1555200Hi,i,i,o,o -77760Hi,i,i,o,i) - 16(1 +9^)Ho,-i,o^2 

+ y(l-2x)Ho,o,-i,o,o + ^(l-^)(12Hi^2V84Hi,o^3 + 24Hi,i^3-96Hi,o,-i^2 
+ 72Hi,o,o^2 + 24Hi, 1,0^2 — 7Hi, 1,1,0 + 20Ho,-i, -1,-1,0 — 14Ho,-i,o,-i,o — 2Ho,-i,o,i,o 

- 26Ho,o,- 1,-1,0 + 2Ho,o,i,o,i — 2Ho,o,i,i,o + 48Hi,o,-i,o,o + 96Hi,o,-i,o,i + 24Hi,o,o,-i,o 
-24Hi,o,o,o,o -48Hi,o,o,o,i +48Hi,i,o,-i,o -24Hi,i,o,o,o) - ^(1 +^)(12H-i^2^ 
-720H-i,-i^3 + 360H-i,o^3+960H-i,-i,-i^2-600H-i,-i,o^2-480H-i,o,-i^2 

+ 420H-i,o,o^2 - 50Ho,i,i^2 + 720H-i,o,i,o + lOOOH-1,0,1,1 + 10Hi,o,i,o -480H-i,-i,-i,o,o 

- 960H- 1 ,- 1 ,- 1 ,0,1 + 120H- 1,-1 ,0,0,0 + 480H- 1,-1 ,0,0, 1 + 240H- 1 ,0,-1,0,0 + 480H- 1 ,0,- 1 ,0, 1 
+ 120H- 1,0,0,- 1,0 — 120H- 1,0,0,0,0 — 300H-i,o,o,o,i + 60H- 1,0,1, 0,0 +45Ho,o,o, 1,0 

+ 75Ho,o,o,i,i + 5Ho,o,i,o,o — 75Ho,i,o,o,i — 5Ho,i,o,i,o + 105Ho,i,i,o,o) +7jc)Ho,i,o,o,o 

8 16 2 4 

- -(2 + 5^)H-i,o,-i,o - y (1 + 13^)Ho,i,o,-i,o - -(1 +33^)Ho,i^2 + -(1 + 35^)Ho,o,-i^2 

- 8(6 - 5;c)Ho,- 1,-1,0 + 4(41 - 8x^1 - 39x + 8x2 )Hi, 0,0,0 

4 8 8 

- -(2 + 3x)Hi,o,i,i + 3(1- 5x)Ho,o,o,-i,o - 3 (3 - 24x-i - 5x + 24x2)Hi,o,-i,o 

4 4 

+ - (4 - 27x- ' - 14x + 27x2 )Hi ,0^2 _ _ (5 _ 5 1 _ _^ 5 l;<:2)Hi,o,o, 1 

16 8 4 

+ _(5 + 43x)Ho,-i,-i^2-3(5 + 91x)Ho,-i,-i,o,i+8(4 + 3x)Ho,o,o^2+3(5-llx)Hi,i^2 



30 



8 8 8 32 

+ -(7 + 5;c)Ho,_i,o,o,o " 3 + 1 l^)Ho,o,-i,o,i - 3 (5 + 17;c)Ho,i,o^2 + y (5 + 2jc)Ho,i^3 

+ ^ (7 + 89^)Ho,- 1 ,0,0, 1 + ^ (9 - 4^)Hi, 1 ,0, 1 - ^ ( 1 1 + 9^)Ho,o,o,o,i +^(11 + 15^)Ho,o, 1 ^ 2 

+ ^(14+ 15a:)H_i,_l_i,o + y (17 + 12;^-' + 17;c+ 12;c2)H_i,_i,o,i + ^(17 + 180;c)Ho,o,i,i 

4 2 2 

- -(19+ 125;c)Ho,_i^3 + ^(21 - 1315;c+ 128jc2)Ho^2 - ^(4363 - 5244;c + 2742;c2)Ho,o 

- ^ (29 + 67;c)Ho,- 1,-1,0,0 + ^(30 + 264;^"^ + 572;c - 67;c2)Ho,-i,o + 2(7 + 73;c)Ho,o,i,o 

- 16(31 + Sx^^ +3lx + 8;c^)H_i,o,o,i + j(25 - 5 - 29;c + 5 1a;^)Hi, 1,0,0 
+ ^(34 + 75^;^^ -20;c-75;c2)Hi^3-^(56- 197;c"^ +65;c- 197;c2)H_i,o,o 

+ ^(39 + 1315^ - 262^2 ^Ho,o,i + ^(34 - 37^)Hi,o,o + ^(37 + 63^)Ho, 1,0,1 " ^(49 + 15^)^5 
-4(49 + 8;c-^ +47;c + 8;c2)H_ 1,0,0,0 - -^(51 + 20;c)Ho^2^ - ^(67 + 27;c + 48;c2)Ho,o,-i,o 

- ^(76 - 12481;c) - ^(77 + 38^:+ 198;c2)Hi,o,i + ^(77 + 41^:+ 198;c2)Hi,i,o 
ol 9 9 

- ^(83 + 5^)^2^3 + ^(84- 16^ + 51^2 J ^^^^ _ ^(131 +9^ + 48^2)jj^ 

8 2 1 

- -(224 -9x + 96x2)Ho,_i,o,i + -(237 + 271x + 396x2)Ho,i,o + -(2533- 1915a;)Hi,i 

- ^(258 + 192x-' +257X+ 192x2)H_i,_i^2 + ^(302 + 144x-' +301x+ 144x2)H_i^3 
+ ^ (3097 + 1 3 17x- ' + 3 15 Ix + 1 37 Ix^ )H_ 1 ,0 - ^ (348 - 67x- ^ + 3 1 8x - 67x2 )H_ ^ 

- ^(788 + 67x"' +788x + 67x2)H_i,o,i + ^(373 + 120x"' +372x+ 120x2)H_i,o^2 

- ^(408 + 482X+ 197x2 )Ho,o,o + ^(430 - 3x+ 192;c2)Ho,_i^2 + ^(977 - 1223x)Hi,o 
_ ^ (449 + 40 Ix + 1 86x2 ^ jj^^^^^ + ^ (3 15 + 49 u + 90x2 ^ jj^^^^^ 1-^(31 + 63x)Ho, 1 ,1,0 

- ^(721 +2035X+ 1523x2)^3 _ ^(1095 - 619x+ 1296x2)^2^ + ^(20601 - 5183x)Ho 

9 15 162 

+ ^(1228 + 201x-^ + 1258x + 201x2)H_i^2 + ^(502 + 67x-^ - 242x + 329x2)Hi^2 

- ^(367 + 5x + 219x2)Ho^3 + ^(3459 - 2365x + 5484x2)^2 + ^(8929 + 2365x)Ho,i 

1 8 1 1 f 1909753 

+ —(40039 - 23821x)Hi + -(19Ho,o^3 - 6Ho,o,oAo) + t^i (x) + -g2{x) + 



416„ 105739„ 248„ 2 143282„ ^ 25184„ 2 3512„ 3K. , 

+ CpC^nf (2700^3 + 830^2 + 1026^2^ - 1440H_i^3 - 3000H_i^2 
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- 1305Ho + 990HoC3 + 345Ho^2 + 1800H_i,_i^2 + 1660H_i,o - 1260H_i,o^2 

- 2340Ho -1^2 - 1910Ho,o + 450Ho,o^2 - 3600H_i _ 1,0 + 4800H_ 1,0,0 + 1200H_i,o,i 
+ 3000Ho,-i,o - 1245Ho,o,o - 600Ho,o,i + 3600H_i,_i,_i,o - 6120H_i,_i,o,o 

- 3960H_i, 0,-1,0 + 6120H_i,o,o,o + 360H_i,o,o,i - 720H_i,o,i,i - 3960Ho,- 1,-1,0 

+ 6120Ho,_i,o,o + 360Ho,-i,o,i +4500Ho,o,-i,o -2070Ho,oao -360Ho,oai +360Ho,o,i,i) 

+ ^^pqq (a:) (402265 + 66690^3 - 521010^2 + 27216^2^ + 895905Ho + 1 19880Ho^3 

- 210870Ho^2 + 338895Hi + 12960Hi ^3 - 168480Hi ^2 - 77760Ho,_i ^2 + 731250Ho,o 

- 74520Ho,o^2 + 491 130Ho,i - 103680Ho,i^2 + 193320Hi,o - 139320Hi,o^2 + 209520Hi,i 

- 81000Hi,i^2 + 126360Ho,_i,o + 318060Ho,o,o +220050Ho,o,i + 174150Ho,i,o 
+ 142560Ho,i,i + 119880Hi,o,o + 103680Hi,o,i +38880Hi,i,o + 71280Hi,i,i 

+ 155520Ho,_ 1, 0,0 + 77760H0,- 1,0,1 + 155520Ho,o,-i,o + 74520Ho,oao +71280Ho,o,o,i 
+ 29160Ho,o,i,o + 35640Ho,o,i,i + 22680Ho,i,o,o + 32400Ho,i,o,i - 32400Ho,i,i,o 
+ 77760Hi,o,-i,o + 51840Hi,o,o,o + 106920Hi,o,o,i -29160Hi,o,i,o +29160Hi,o,i,i 
-42120Hi,i,o,o + 16200Hi,i,o,i -45360Hi,i,i,o) + ^(1 -^)(5Hi,i^2 - 5Hi,o,o -4Ho,-i,o,i 

- 4Hi,o -1,0 - 6Hi,o,o,o + Hi,o,i,o - Hi,i,o,i + Hi,i,i,o) - - (5 + Ix^ -5x-2x )Hi,o,o,i 

-^(l+^)(12H_i^3-12H_i,_i^2 + 6H_i,oC2-7Ho,i,i-24H_i,_i,_i,o + 36H_i,_i,o,o 

8 16 
+ 24H_ 1 ,0,- 1 ,0 - 1 8H_ 1 ,0,0,0 + 24Ho,- 1 1 ,0) + — ( 1 + 4^ + 8^^ ) ^2 V y ( 1 + 7^)Ho,- 1 ,0,0 

-y(2-A;"^+2x-x2)H_i,_i,o + ^(3-2x"' -3a; + 2,x2)Hi^3-^(3-35;c + 72,x2)Ho,i,o 

+ I (4 - 5^)Ho, 1 C2 + I (4 + 3^ + 2x2)Ho^3 + ^ (5447 - 1 372 - ( 1 - ^ + jj^^ ^ ^^^^ 
5 J ol 3 

+ ^(7 + 2x-i - 7;c - 2x2)Hi,o^2 - y (8 + 3;c-^ + 8;c + 3;c2)H_i,o,i + ^(20 - 77;c)Hi,i 

32 16 4 

- y (8 + 9x-i + 8;c + 9;c2)H_i 0,0 + y (10 + 9x-^ + lOx + 9x2)H_i^2 + ^ (37 - 61;c)Hi,o 

4 2 2 

- - ( 1 1 - 1 2x" ' + 20x + 48^:^ )H 1 ^2 + ^ ( 1 69 - 1 1 3x + 96x2 ) ^3 + — (667 - 67 lx)Ho, 1 

+ l(35_4x + 36x2)(Hi,o,i-Hi,i,o) + ^(81 + 52x + 36x2)Ho,o,o + y(^"^-^^)Hi,i,o,o 

8 8 

- 27(158 + 141x"^ + 158x+ 141x2)H_i,o - -(15 + 24x"' - llx+ 12x2)Ho,_i,o 

2 2 1 

+ -(207+ 149x+ 120x2 )Ho,o,i - ^(267 + 149x+ 168x2 )Ho^2 + — (7651 - 23577x)Ho 

- ^(1299 -671x + 564x2)^2 + ^(1685 -1219x + 564x2)Ho,o + ^(1855 -4693x)Hi 

8 8 

~ o (^ + 2x2)(Ho,o^2 — Ho,o,o,i) — o (8Ho,-i^2 — 12Ho,o,-i,o + Ho,o,i,o — Ho,i,o,i +Ho,i,i,o)x 
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18314 
81 
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Eq. (IA.4I ) includes, with the colour factor C^(C4 — 2 C^)^, the sum of the inverse Mellin transforms 
of the functions gj (N) and g2{N) in Eqs. (12.221) . These functions are given by 

gi (x) = -2(1 - (1 -x)-i)(4H_i^2 - 2H_i,o,o - 4H_i,o,i + Ho,o,o + 2Ho,oj - 3HoC2 - 3^3) 
+ 2((1 - _ (1 (4Ho,-i C2-2Ho,-i, 0,0 -4Ho,- 1,0,1 +2Ho,o,o,i - 3Ho,oC2 

-3HoC3 + Ho,o,o,o-^C2')+2(l-(l+x)"^)(2Ho,-i,o-Ho,o,o + HoC2 + 2C3) 

-((l+x)-2_(i+^)-i)(^4Ho,o,-i,o + 2Ho,oC2+4HoC3-2Ho,o,o,o + yCl) , (A.5) 

g2{x) = -2(2(1 -3(1 -x)-2 + (l-x)-i)(4Ho,-iC2+Ho,o,o,o-2Ho,-i,o,o 

-4Ho,-i,o,i +2Ho,o,o,i - 3Ho,oC2 -3HoC3 - - 2^2(1 - (1 -x)-') 

_4((l_^)-2_(i_^)-i)(4H_i^2-2H_i,o,o-4H_i,o,i+2Ho,o,i+Ho,o,o-3HoC2-3C3) 

+ (2(l+.^)-3_3(l+x)-2 + (i+^)-i)(^4Ho,o,-i,o-2Ho,o,o,o + 2Ho,oC2 + 4HoC3 + yCl) 
+ 4 ( ( 1 + x) -2 _ ( 1 + x) - 1 ) (2Ho,- 1 ,0 - Ho,o,o + H0C2 + 2^3 ) 

-2(l-(l+x)-i)(2H_i,o + 2Ho,i-C2) + (C2-C3)5(l-x) . (A.6) 

Note that the 1 / ( 1 — x) terms in Eqs. (IA.5I) and (IA.6I) do not imply the presence of +-distributions, 
but large-x cancellations between the harmonic poly logarithms in the round brackets. These can- 
cellations become numerically problematic for I — x <^ I, a region unavoidable in moment and 
Mellin convolution integrals, even with the high-accuracy code of Ref. [47]. In this region one can 
instead use 

Ui+^2)W = (C2-C3)5(l-x)-(l-x){ln(l-x)-F|-C3-^} 

-^-{l-xf + 0[{l-x)']. (A.7) 



Appendix B: Subleading large-x contributions to F2 



ns 



In this final appendix we consider the subleading (integrable) large-x logarithms ln*^(l — x) for 
the quark contributions to F2 and their relation to the coefficients for F3 given in Eqs. (13.91) - 
(13.161) above. The corresponding three-loop results for (where these terms are the leading 
contributions) have been given in Eq. (4.31) - (4.34) of Ref. [18]. At two loops we have 
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The corresponding contributions to the third-order coefficient function are given by 
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A comparison of the large-x limit of the complete charged-current quark coefficient functions 
Ci,2,3,L, defined analogous to Eq. (13.11) above, now reveals that 

Ci(x,as) = C2{x,as)-CL{x,a,) = C3{x,a,) + o{l -x) (B.9) 



to the third order in tts. I.e., not only are the above log-enhanced terms of Ci and C3 related, 
but also the constants for x ^ I, leaving the difference vanish as 1 /N^ in moment space. To the 
best of our knowledge, Eq. (IB. 91 ) has not been noted so far in the literature even on the level of 
the two-loop coefficient functions of Refs. [4-6]. In this context it appears worthwhile to recall 
that the non- singlet coefficient function for the polarized structure function [no relation to the 
quantity in Eq. (IA.5I) 1 is identical to C3 up to the d"^^'''dabc term irrelevant at large x, and that an 
analogous 1 /N^ suppression of the helicity flip for x ^ 1 is present in polarized splitting functions, 
see Ref. [61] and references therein. 
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